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light.  DNA  was  not  damaged  at  any  energy  density  available  including  levels 
sufficient  to  cause  dielectric  breakdown  in  water.  Experimental  studies  of 
electromagnetic  stress-induced  bi refri ngence  in  DNA  base-pairs  were  also  carried 
out  in  an  attempt  to  establish  a  lower  limit  on  the  restraining  forces  governing 
tilting  of  the  DNA  bases  with  respect  to  the  helix  axis.  The  irradiation 
experiments  at  1060  nm  with  the  bilipid  layer  vesicles  indicated  a  damage 
threshold  of  ~600  mJ/cmf  for  entire  pulse  trains  of  -100  pulses  and  ~9  inJ/cm2 
for  single  pulses*  as  determined  by  dynamic  light  scattering.  In  both  the 
multiple-  and  single-pulse  cases,  the  peak  optical  electric  field  incident  c-n 
the  vesicles  was  6-7xl05  V/cm,  about  an  order  of  magnitude  above  the  membrane 
potentials.  -It  appears  likely  that  direct  electrostricti ve  forces  disrupted  the 
vesicle  membranes  and  facilitated  their  transition  to  stacked  lamellar  struc¬ 
tures  known  as  liposomes.- 

Retinal  damage  thresholds  in  the  Maeaca /fascicular is  were  determined  for 
irradiation  with  single  ultrashort  pulses  and  entire  pulse  trains  of  2nd  harmon¬ 
ic  light.  In  terms  of  irradiance  at  the  retina,  the  fundoscopically  deter¬ 
mined  24-hr  postexposure  thresholds  were  4.4  mJ/cm2  and  540  mJ/cm2  respectively. 
The  peak  electric  fields  in  both  cases  were  of  the  order  of  6-7xl05  V/cm2,  as  in 
the  case  of  the  vesicles.  Disruption  of  the  cellular  membranes  is  suggested  as 
the  threshold  damage  mechanism.  In  addition,  it  is  postulated  that  for  pulse 
trains,  irreversible  damage  occurs  very  early  in  the  train,  at  the  first  pulse 
to  attain  or  exceed  the  threshold  electric  field. 

The  corneas  of  the  same  primates  used  in  the  retinal  studies  were  irradiated 
with  4th  harmonic  (265  nm)  mode-locked  pulse  trains  derived  from  the  NdrGlass 
laser.  Fluorescein  slit-lamp  examinations  at  24-hr  post  exposure  revealed  a 
damage  threshold  of  8.2  mJ/cm2.  All  damage  was  confined  to  the  corneal  epi¬ 
thelium.  It  is  postulated  that  photochemical  processes  such  as  coagulation  or 
denaturation  of  nucl eoproteins  and  nucleic  acids  govern  the  damage  mechanism 
in  this  case. 
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BIOLOGICAL  DAMAGE  THRESHOLD  INDUCED  BY  ULTRASHORT  FUNDAMENTAL,  2ND, 
AND  4TH  HARMONIC  LIGHT  PULSES  FROM  A  MODE-LOCKED  Nd: GLASS  LASER 


INTRODUCTION 

The  effects  of  laser  radiation  on  various  biological  systems  have  been 
studied  extensively  for  several  years  (1-5).  The  eye  has  been  found  to  be 
particularly  susceptible  to  damage  by  exposure  to  radiation  ranging  from  the 
ultraviolet  (UV)  to  the  infrared  (IR)  portions  of  the  optical  spectrum  (1-3, 
6-31).  Depending  on  the  wavelength  of  the  laser  radiation,  ocular  damage  is 
most  likely  to  occur  in  the  cornea,  the  lens,  or  the  chorioretinal  tissues  of 
the  eye.  The  clear  tissues  of  the  eye  are  quite  transparent  in  the  visible 
and  near  IR  (32);  consequently,  radiation  in  this  wavelength  range  affects 
primarily  the  retina.  In  the  UV  the  bulk  of  the  incident  radiation  is  absorbed 
in  the  anterior  portion  of  the  eye;  below  about  300  nm,  the  absorption 
occurs  entirely  in  the  cornea,  usually  within  the  corneal  epithelium  (2,25). 

It  follows  then  that  ocular  damage  from  the  visible  and  near  IR  is  sustained 
primarily  in  the  retinal  tissues;  while  from  the  UV,  it  is  sustained  primarily 
in  the  corneal  tissues. 

The  mechanisms  responsible  for  UV-induced  damage  in  the  cornea  are 
believed  to  be  predominantly  photochemical  (25,27).  However,  relatively  little 
work  has  been  done  on  the  ocular  hazards  of  UV  laser  radiation,  particularly 
in  the  ultrashort  pulse  regime  and  with  respect  to  the  effects  of  pulse  dura¬ 
tion  in  general. 

By  contrast,  a  large  body  of  data  exists  on  the  effects  of  visible  and 

near-IR  laser  radiation  on  the  retina  (6-24,  28-31).  The  bulk  of  the  work 

has  been  carried  out  with  CW  or  pulsed  lasers  of  various  wavelengths  in  the 

nanosecond  or  longer  pulse  duration  regime.  Until  recently,  relatively  little 

-12 

attention  was  paid  to  the  ocular  hazards  of  ultrashort  laser  pulses  (10 

-11  -8 
10  sec)  (28-31).  For  the  longer  pulse  durations  (>10  sec),  the  injury  to 

the  retinal  tissue  is  believed  to  be  caused  by  local  temperature  rise  and 

resulting  protein  denaturation  and  enzyme  inactivation.  Typical  threshold 

values  of  radiant  exposure  at  the  retina  that  cause  observable  lesions  approach 
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2  -8 
^1  J/cm  as  the  exposure  time  is  reduced  to  10  sec.  In  the  picosecond  re¬ 
gime,  there  is  less  agreement  on  the  damage  thresholds,  and  values  as  low  as 
2xl0~3  J/cm^  (30)  and  as  high  as  2  J/cm^  (28)  have  been  quoted. 

There  is  general  agreement,  however,  that  in  the  picosecond  regime  the 
tissue  damage  is  not  of  thermal  origin  (4,28-31).  Several  possible  damage 
mechanisms  have  been  suggested,  such  as  intense  acoustic  transients,  direct 
breakdown  in  the  bulk,  multiphoton  ionization,  and  free  radical  formation 
(4,24,28-31).  Other  possible  mechanisms  include  direct  electromagnetic  field- 
induced  stresses  in  the  nucleic  acids,  proteins,  and  cell  and  vesicle  lipid 
membranes  (30,31,33,34). 

Identification  of  specific  damage  mechanisms  in  living  tissue  is  diffi¬ 
cult  because,  unlike  the  optically  pure  samples  used  in  inanimate  material 
studies,  biological  materials  are  heavily  concentrated  with  optical  inclusions, 
free-charge  regions,  and  dielectric  absorption  discontinuities.  Furthermore, 
the  end  point  resulting  from  energy  deposition  may  be  reached  by  a  complicated 
set  of  pathological  changes,  with  a  result  that  is  difficult  to  interpret. 

The  important  constituents  of  a  living  cell  that  are  likely  candidates  for 
sites  of  radiation  damage  are  its  proteins,  including  fibrous  structural  and 
contractile  proteins,  enzymes,  and  other  polypeptides;  its  nucleic  acids, 
which  are  generally  found  in  cells  to  be  complexed  with  basic  proteins  (his¬ 
tones)  in  a  form  called  chromatin,  or  complexed  with  polyamines;  and,  finally, 
its  cell  and  vesicle  lipid  membranes  (34). 

We  have  approached  this  problem  by  (1)  investigating  the  effects  of  ultra¬ 
short-pulse  laser  radiation  on  isolated  biological  macromolecular  systems 
similar  to  constituents  of  living  cells  that  are  susceptible  to  damage  and  (2) 
comparing  the  results  to  data  obtained  from  direct  ocular  damage  experiments. 

The  first  phase  of  the  program  described  in  this  report  was  designed  to 
study  the  effect  of  ultrashort  pulses  of  1060-nm  light,  and  its  2nd  harmonic, 
derived  from  a  mode-locked  Nd:Glass  laser,  on  aqueous  suspensions  of  macromo¬ 
lecular  structures  such  as  calf-thymus  DNA,  poly(L-lysine) ,  DPPC  vesicles,  and 
egg-yolk  lecithin  vesicles.  The  experiments  have  encompassed  a  wide  range  of 
pulse  energies  for  both  single  pulses  and  entire  trains  of  mode-locked  pulses. 
Particle  damage  has  been  monitored  by  one  or  more  of  the  following  techniques 
(depending  on  the  particular  macromolecule):  dynamic  light  scattering,  gel 
electrophoresis,  and  low-shear  viscometry.  In  the  case  of  DNA,  which  suffered 


no  damage  detectable  by  these  techniques  at  the  energy  densities  attainable, 
we  have  also  carried  oat  a  theoretical  and  experimental  investigation  of 
transient  distortion  birefringence  to  ascertain  if  any  detectable  internal 
strain  is  induced  by  the  laser  pulses. 

The  second  objective  of  our  program  was  to  determine  direct  ocular  dam¬ 
age  thresholds  in  primates  resulting  from  the  2nd  harmonic  (530  nm)  and  the 
4th  harmonic  (265  nm)  ultrashort-pulse  radiation  derived  from  the  same  laser, 
and  to  compare  the  results  with  those  of  the  macromolecular  studies,  with  a 
view  to  identify  more  closely  the  most  likely  damage  mechanism(s) .  These 
studies  were  also  aimed  to  help  establish  new  laser  safety  standards  in  the 
picosecond  time  regime  at  these  wavelengths. 

This  report  is  divided  into  two  major  sections:  Part  I  deals  with  the 
macromolecular  irradiation  studies,  and  Part  II  treats  the  experimental  de¬ 
termination  of  ocular  damage  thresholds  in  primate  eyes. 
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PART  I:  MACROMOLECULAR  DAMAGE  STUDIES 


REVIEW  OF  PRIOR  WORK 

Important  macromolecular  constituents  of  mammalian  cells  that  are  likely 
candidates  for  initial  sites  of  ultrashort-pulse  laser  radiation  damage  are: 

(a)  their  proteins,  including  fibrous  structural  and  contractile  proteins, 
enzymes,  and  other  polypeptides;  (b)  their  nucleic  acids,  which  are  general¬ 
ly  found  in  cells  to  be  complexed  with  basic  proteins  (histones)  in  a  form 
called  chromatin,  or  complexed  with  polyamines  such  as  spermidine;  and  (c) 
their  noncovalent  lipid  membrane  structures,  including  those  of  the  outer 
integument,  various  vesicles,  mitochondria,  and  other  endoplasmic  inclusions. 

Work  described  in  our  previous  technical  report  on  this  subject  (34)  in¬ 
dicated  that  the  polypeptide  poly(L-lysine)  was  unaffected  by  either  single 

o 

ultrashort  pulses  of  1060-nm  laser  light  with  energies  up  to  1.37  mJ/cm  or 

successions  of  up  to  four  pulse  trains,  each  consisting  of  '.100  pulses  and 

o~ 

having  integrated  energy  densities  of  up  to  150  mJ/cm  .  It  was  inferred  that 
polypeptides,  or  proteins,  are  unlikely  to  be  the  initial  sites  of  any  bio¬ 
logical  damage  at  these  or  lower  energy  densities. 

The  major  effort  of  our  previous  investigations  was  to  study  the  effects 
of  ultrashcrt  pulses  of  1060-nm  laser  radiation  on  DNA.  Three  potential  dam¬ 
age  mechanisms  were  theoretically  analyzed  for  their  resulting  distributions 
of  fragment  sizes.  We  concluded  that  a  clear-cut  distinction  between  damage 
mechanisms  could  in  principle  be  made  in  the  event  that  damage  was  unequivo¬ 
cal  ly  observed.  Unfortunately  for  us,  but  perhaps  fortunately  for  laser-damage 
victims,  breakage  of  the  DNA  was  unambiguously  observed  in  only  one  sample,  a 
result  that  we  now  strongly  suspect  was  spurious.  Most  of  the  experiments 
revealed  no  detectable  damage.  Serious  problems  with  both  radiated  and  con¬ 
trol  DNA  samples  toward  the  end  of  the  contract  period  resulted  in  a  number 
of  ambiguous  experiments  in  which  the  DNA  exhibited  an  anomalously  high  affin¬ 
ity  for  the  Mi  1 1 i pore  filters,  with  consequent  slow  filtration  and  loss  of 
sample  and  subsequent  loss  of  precision  in  the  dynamic  light  scattering 
measurements,  which  gave  no  positive  indication  of  damage  in  any  case.  Because 
of  these  technical  difficulties  arm  because  sf  thr  importance  of  observing 


directly  the  fragment  distribution,  we  undertook  the  development  of  a  new 
analysis  technique,  gel -electrophoresi s .  However,  final  completion  of  the 
necessary  apparatus  and  application  of  that  technique  to  study  irradiated  DNA 
samples  had  to  await  the  present  contract  period. 

In  summary,  no  detectable  damage  of  DNA  was  produced  by  single  pulses  up 

9  ? 

to  0.75  mJ/cm'",  or  two  pulse  trains  up  to  230  mJ/cm  total  eneray  density,  or 

2 

four  pulse  trains  up  to  440  mJ/cm  .  At  higher  total  energy  densities,  the 
results  were  ambiguous  due  to  technical  problems.  A  primary  objective  of  the 
present  work  has  been  to  obtain  unambiguous  results  for  DNA  irradiated  at  the 
highest  obtainable  total  energy  densities. 

Some  initial  work  was  performed  on  the  physical  characterization  of  a 
preparation  of  purple  membranes  from  Hal obacteri urn  halobium,  which  we  hoped 
would  provide  a  tractable  model  for  a  retinal  membrane-protein  system.  Un¬ 
fortunately,  the  ul trasonically  dispersed  suspensions  of  this  material  were 
neither  as  reproducible  nor  as  stable  as  desired,  so  no  irradiation  experi¬ 
ments  were  undertaken. 


OBJECTIVES 

At  the  outset  of  the  present  investigations  we  intended  to  study  DNA  ir¬ 
radiated  at  the  highest  attainable  energy  densities--using  dynamic  light 
scattering,  low-shear  viscometry,  and  especially,  gel  electrophoresis--to 
detect  any  damage  sustained  and  to  analyze  the  distribution  of  product  frag¬ 
ments.  When  it  became  clear  that  DNA  was  not  detectably  damaged  by  even  the 
most  energetic  available  laser  pulses  or  whole  trains  of  such  pulses,  we  de¬ 
cided  to  investigate  the  transient  distortion  birefringence  induced  by  such 
picosecond  pulses.  Our  idea  was  that,  even  though  rupture  of  covalent  bonds 
was  not  observed,  we  might  still  be  able  to  observe  transient  distortion  in¬ 
duced  by  the  high-power  optical  pulse;  however,  no  transient  distortion  bire¬ 
fringence  was  ever  observed. 

When  it  became  clear  that  the  covalent  bonds  of  both  polypeptides  and 
nucleic  acids  were  unaffected  by  ultrashort  pulses  of  1060-nm  light  at  the 
energy  densities  attained  in  these  studies,  we  decided  to  focus  our  effort  on 
noncovalent  lipid  membrane  vesicles.  Indeed,  it  was  unambiguously  and  repro- 
ducibly  observed  that  irreversible  changes  were  induced  in  suspensions  of  di- 

palmitoyl  phosphatidyl  choline  vesicles  by  trains  of  pulses  with  total  enemy 

2 

densities  exceeding  600  mJ/cm  . 
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SEARCH  FOR  PICOSECOND  OPTICAL  STRESS-INDUCED  DAMAGE  IN  DNA 


The  quality  control  problems  with  the  calf  thymus  DNA  and/or  filtration 
system  that  surfaced  in  the  final  phase  of  the  early  work  have  been  essenti¬ 
ally  overcome  by  a  combination  of  procedures,  includinq 

(1)  washing  the  filters  in  a  hot  solution  containing  0.0115  M  Na2C03  and 
0.0015  EDTA,  and 

(2)  working  with  the  DNA  in  a  higher  pH  buffer,  1  M  NaCl,  0.01  M  EDTA, 
and  0.05  M  Na^CO^,  pH  9.3. 

We  do  not  understand  why  DNA  passes  through  Millipore  filters  so  much  more 
easily  at  high  pH  (>  8.5)  than  at  neutral  pH,  but  that  has  been  a  repeated 
observation  in  this  laboratory. 

Laser  experiments  with  DNA  at  1060  nm  were  performed  following  the  pro¬ 
tocol  described  in  our  previous  report  (34)  with  a  somewhat  modified  laser 
system.  Descriptions  of  the  irradiation  apparatus  and  experimental  protocol 
are  reproduced  here  in  Appendixes  A  and  8  respectively.  We  temporarily  re¬ 
moved  the  optics  used  for  pulse  chronography  in  order  to  increase  the  pulse 
energy  available  at  the  irradiation  sample  cell.  The  laser  was  operated 
multimode  with  the  larger  diameter  laser  beam  focused  down  with  a  lens  system 
to  match  the  ID  of  the  sample  cell.  This  further  increased  the  irradiation 
intensity.  The  beam  profile,  although  not  Gaussian,  was  free  of  large-scale 
intensity  fluctuations  (i.e.,  hot  spots).  Spatial  intensity  ripples  were 
within  10"  of  the  mean  local  intensity. 

Instrumentation  and  Techniques  for  Studying  Possible  Damage  to  DNA 

Gel  Electrophoresis--A  gel  electrophoresis  apparatus  was  constructed 
using  a  J-wick  flat-bed  configuration  designed  for  very  dilute  (0.5°'  or  less) 
agarose  gels.  The  entire  unit  was  housed  in  a  Lucite  safety  box  with  micro- 
switches  in  the  latch  to  guarantee  interruption  of  the  high-voltaqe  leads 
whenever  the  box  was  open.  Considerable  experimentation  was  done  with  methods 
of  preparing,  pouring,  storing,  loading,  running,  and  developing  the  dilute 
agarose  gels.  Some  of  the  most  important  guidelines  to  emerge  from  this 
activity  were  the  following.- 


(1)  Agarose  concentrations  of  at  least  0.3%  (w/V)  were  required  to  form 
gels  strong  enough  to  permit  the  necessary  manipulation  without 
breaking. 

(2)  A  "frame"  several  mm  wide  and  1-2  mm  deep  of  more  concentrated  gel 
(1%)  on  the  periphery  was  essential  to  provide  adequate  mechanical 
stability  for  manipulation  of  the  more  dilute  bulk  gel. 

(3)  Interchanging  the  buffer  in  the  anode  and  cathode  compartments 
every  few  hours  was  effective  in  moderating  pH  drift. 

(4)  Forming  the  gel  on  top  of  a  removable  plastic  slab  greatly  facili¬ 
tated  the  subsequent  manipulations. 

Typically,  10-20  microliters  of  solution  containing  0.03-0.10  mg/ml  DNA 
were  loaded  into  one  of  the  wells,  or  slots,  on  the  "starting  line"  of  a  cold 
gel.  These  wells  were  impressions  of  the  embedded  teeth  of  a  suspended  comb 
remaining  after  t.ie  gel  had  cooled.  The  loaded  wells  were  capped  with  a  drop 
of  warm  agarose  solution  that  was  allowed  to  harden.  The  electrophoresis 
buffer  that  we  have  found  works  best,  and  gives  the  least  pH  drift,  is  40  mM 
Tris  (base),  1  mM  NaEDTA,  and  5  mM  H3B03.  The  voltage  was  initially  set  at 
60  kV  for  20-30  min  until  the  DNA  penetrated  into  the  gel  (from  the  well), 
then  was  maintained  at  30  kV  thereafter,  which  produced  about  2  mA  current. 
Bromthymol  blue  was  added  to  one  of  the  wells  as  a  migration  indicator.  It 
moves  with  two  to  three  times  the  speed  of  the  native  intact  DNA.  Typical 
run  times  were  18-24  hr.  The  gel  slab  was  then  removed  from  the  apparatus  and 
soaked  in  a  solution  containing  ethidium  bromide  (%0.001  mg/ml)  for  20-30  min. 
Then  the  gel  was  transferred  to  the  illuminating  upper  surface  of  a  UV  black- 
light  box.  Ethidium  cation  fluoresces  with  a  quantum  efficiency  about  200  x 
higher  when  bound  to  DNA  than  when  free  in  solution,  thus  its  bright  orange 
fluorescence  locates  the  position  of  the  DNA  "bands"  in  the  gel.  These  bands 
could  be  photographed  using  an  appropriate  filter  to  block  the  exciting  UV 
light.  The  resulting  negatives  could  be  scanned  by  densitometer  to  provide 
a  very  precise  quantitative  picture  of  the  DNA  distributions  in  the  bands. 
Although  this  was  done  for  a  few  runs,  using  equipment  in  the  laboratory  of 
W.  Fangman,  direct  visual  inspection  generally  sufficed  to  demonstrate  that 
no  damage  was  occurring.  Therefore,  densitometry  was  not  further  pursued, 
and  that  subject  will  not  be  discussed  in  the  sequel. 

At  the  time  of  writing,  the  total  number  of  DNA  samples  analyzed  by  gel 
electrophoresis  in  this  laboratory  (for  this  and  other  projects)  is  well  over 
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100.  We  have  found  that  the  technique  is  extraordinarily  sensitive  for  de¬ 
tecting  double-strand  breaks,  and  practically,  if  net  totally,  insensitive  to 
single-strand  breaks.  A  surprising  finding  of  our  ongoing  investigations  of 
DNA  was  that  even  the  internal  motions  manifested  in  the  dynamic  light  scat¬ 
tering  were  completely  insensitive  to  single-strand  breaks  at  neutral  pH  (35). 
In  fact,  data  in  the  older  literature  (36)  indicate  that  the  sedimentation 
coefficient  and  radius  of  gyration  obtained  from  static  light  scattering,  both 
of  which  monitor  the  bending  rigidity,  are  insensitive  to  single-strand  breaks 
at  neutral  pH.  Thus,  the  insensitivity  of  the  gel  electrophoresis  to  single¬ 
strand  breaks  is  consistent  with  these  other  observations.  A  recent  discovery 
has  shown  that  single-strand  breaks,  bound  protein  contaminants,  and  bound 
spermidine,  can  all  induce  profound  changes  in  the  Rouse-Zimm  model  internal 
motion  parameters  at  pH  10.2  (37).  These  facts  may  be  useful  in  any  future 
DNA  damage  studies. 

The  presence  of  double-strand  breaks,  which  we  have  often  observed  in 
samples  contaminated  with  nuclease  activity,  is  manifested  by  lower  molecular 
weight  components  migrating,  or  streaking  well  ahead  of  the  main  band  of  the 
intact  DNA.  A  simple  description  of  the  gel  electrophoresis  technique  as 
applied  to  high  molecular  weight  DNAs  is  given  by  Fangman  (38).  ^ 

Low-Shear  Viscometry--A  Crothers-Zimm-type  rotating  float  viscometer 
with  a  shear  of  0.56  sec'1  has  been  used  to  analyze  irradiated  DNA.  The  low- 
shear  specific  viscosity,  as  defined  by 

Solution  . 

.  - - i 

nsolvent 

is  approximately  proportional  to  ,  where  M  =  molecular  weight,  for 

a  fixed  weight  concentration  of  DNA,  and  thus  provides  a  rather  sensitive 
indication  of  the  average  molecular  weight  in  the  sample  (36).  Considerable 
experimentation  was  required  to  evolve  a  procedure  that  would  yield  repro¬ 
ducible  and  reliable  values  of  the  accuracy  required.  Nutation  of  the  float 
greatly  perturbed  the  measured  viscosities,  and  eventually  was  controlled 
only  by  altering  the  configuration  of  the  rotating  magnets,  careful  monitor¬ 
ing  of  the  sample  volume,  and  scrupulous  attention  to  cleanliness. 

The  intrinsic  viscosity  fY]  =  •  /c  =  60  dl/g  observed  for  the  control 
is  close  to  that  expected  for  DNA  molecules  of  molecular  weight  12x10  ,  which 
is  what  good  samoles  of  calf  thymus  DNA  nearly  always  turn  out  to  be. 


i  ■; 


Dynamic  Light  Scattering- -The  experimental  instrumentation  and  tech¬ 
niques,  including  filtration  procedures,  were  exactly  as  described  in  our  re¬ 
cent  report  (34).  (See  also  Appendix  C.)  This  was  by  far  the  most  difficult 
and  time-consuming  assay  for  damage,  and  was  often  fraught  with  problems  be¬ 
lieved  to  arise  from  dust  inevitably  present  in  the  sample  irradiation  cells, 
which  were  very  difficult  to  clean.  Therefore,  this  technique  assumed  a 
secondary  importance  to  the  gel  electrophoresis  technique  and  was  not  run  on 
all  samples,  especially  those  exhibiting  significant  dust  problems. 

Results  of  Damage  Studies  on  DNA  at  1060  nm 

Table  1  presents  the  results  obtained  for  five  separate  samples  of  the 
same  calf-thymus  DNA,  which  were  exposed  to  picosecond  pulse-envelopes  ranging 
from  a  single  pulse  up  to  four  pulse-trains,  each  containing  the  number  of 
pulses  indicated.  The  maximum  energy  densities  are  significantly  higher  than 
those  reported  previously.  The  differences  in  properties  of  the  irradiated 
and  control  samples  lie  well  within  their  respective  experimental  errors  in 
all  cases.  The  positions  and  widths  of  the  gel  electrophoresis  bands  of  all 
the  irradiated  samples  were  identical  to  those  of  the  control,  as  denoted  by 
the  "neg"  result.  There  is  clearly  no  significant  fragmentation  of  DNA  at 
these  power  levels  under  the  present  conditions. 

The  presence  of  unusually  large  amounts  of  dust  and/or  aggregates  in 
these  samples  precluded  accurate  dynamic  light  scattering  measurements,  but 
the  gel  electrophoresis  was  fortunately  unaffected  by  such  contaminants.  In¬ 
terestingly,  the  occurrence  of  Tyndalls  due  to  dust  and/or  aggregates  dimin¬ 
ished  in  the  order  in  which  the  samples  were  introduced  to  the  cell,  suggest- 
ingstrongly  that  the  dust  was  initially  present  in  the  radiation  cell,  which 
was  very  difficult  to  adequately  clean,  and  was  gradually  washed  out  by  suc¬ 
cessive  samples.  The  control,  which  was  run  last,  was  quite  clean. 

Table  2  presents  results  obtained  for  calf-thymus  DNA  at  still  higher 
energy  densities.  In  fact,  in  sample  5  of  that  Table  the  incident  beam  was 
focused  into  the  sample  cell  with  a  55-mm-focal -length  (f.l.)  lens  to  cause 
dielectric  breakdown  in  the  DNA  solution.  The  focal  spot  size  in  that  case 
was  -55  urn.  The  sample  was  irradiated  with  two  pulse  trains  and  was  reposi¬ 
tioned  after  the  first  train  so  that  breakdown  would  be  produced  at  a  differ¬ 
ent  place  in  the  cell.  The  occurrence  of  breakdown  was  evidenced  by  the 
generation  of  a  greenish  spark  2  mm  long  and  a  sharp  audible  acoustic 
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TABLE  1.  SUMMARY  OF  DNA  IRRADIATION  EXPERIMENTS 


Errors  in  [n]  and  D  are  several  percent  in  each  case. 


transient.  Maximum  electric-field  strengths  in  these  focused  pulse  trains  were 

q 

about  6  x  10  V/m.  Again,  both  dynamic  light  scattering  and  gel  electrophoresis 
indicate  that  there  was  no  significant  fragmentation  of  any  of  these  irradi¬ 
ated  DNA  samples.  In  related  noncontract  research,  we  have  found  that  the 
apparent  translational  diffusion  coefficient  D(o=25°)  is  much  more  sensitive 
to  DNA  fragmentation  by  contaminating  nucleases  than  is  D(e=90°). 

It  may  be  concluded  that  there  is  no  detectable  fragmentation  of  DNA  upon 

irradiation  by  single  picosecond  pulses  of  1060-nm  laser  light  with  energy 

2 

densities  up  to  17  mJ/cm  ,  or  pulse  trains  with  energy  densities  up  to 
o 

2.16  J/cm  ,  or  focused  pulse  trains  with  energy  densities  at  focus  up  to 
4.5  x  10  J/cm  ,  the  latter  being  sufficient  to  produce  local  dielectric 
breakdown. 

Damage  Studies  of  DNA  Irradiated  at  1060  and  530  nm  Simultaneously 

Letokhov  (39)  has  suggested  that  multistage  excitation  of  DNA  in  solution 

to  energy  levels  high  enough  to  produce  dissociation,  or  fragmentation,  could 

9  2\ 

be  achieved  by  simultaneous  irradiation  with  very  intense  (^10  W/cm  )  pico- 

11  _  i  2 

second  (VI 0”  -  10  sec)  pulses  of  one  or  more  frequencies,  either  one  of 
which  is  too  low  to  photodissociate  the  molecule  in  a  single  quantum  event. 

This  viewpoint  ignores  the  effects  of  anharmonicity  in  the  vibrational  mani¬ 
fold  and  assumes  a  10-fold  smaller  relaxation  rate  (at  high  excitation  levels) 
than  the  wel 1 -documented  very  rapid  relaxation  of  energy  among  the  numerous 
(coupled)  nondissociative  modes  at  high  excitation  levels  of  large  molecules 
even  in  the  gas  phase  (40).  The  availability  of  a  great  many  more  coupled 
modes  to  accept  the  energy  in  condensed  phases  will  probably  increase  the 
power  requirements  for  such  multistep  photofragmentation  far  above  the  value 
estimated  by  Letokhov.  Nonetheless,  we  decided  to  investigate  the  effect  of 
irradiation  at  1060  nm  and  530  nm  simultaneously.  Two  530-nm  photons  would 
be  equivalent  to  one  at  265  nm,  close  to  the  maximum  of  the  DNA  ultraviolet 
absorption  band  at  263  nm. 

The  experimental  setup  was  similar  to  that  used  for  the  1060-nm  studies, 
except  that  the  dichroic  beamsplitter  tnat  normally  separates  the  IR  and 
green  components  was  replaced  with  an  uncoated-glass  beamsplitter.  Thus,  the 
2nd  harmonic  could  reach  the  target  cell.  Relative  measurements  of  IR  and 
green  pulse  energy  tor  two  experiments  yielded  an  SHG  conversion  efficiency 
of  3.'  in  one  case  and  IT  in  the  other.  The  DNA  samples  were  contained  in 
the  5-mr-in  x  20-mm  quartz  cells.  Tee  laser  aperture  was  set  at  8  mm  and  a 


50-cm-f . 1 .  lens  was  used  to  reduce  the  beam  diameter  to  5  mm.  By  measuring 
burn  spots  on  Polaroid  film  at  the  location  of  the  target  cell,  we  determined 
that  the  diameter  of  the  2nd  harmonic  beam  was  about  half  that  of  the  IR. 

This  is  due  partly  to  the  square-law  harmonic  generation  and  the  different 
effective  focal  length  of  the  focusing  lens  for  the  two  different  wave¬ 
lengths  . 

Table  3  contains  the  results  of  these  simultaneous  (1060  nm  +  530  nm) 
irradiation  experiments.  Again,  it  may  be  concluded  that  no  fragmentation 
of  the  DNA  was  detected.  Although  sample  1  seems  to  be  giving  slightly  high¬ 
er  apparent  diffubion  coefficients,  its  gel  electrophoresis  band  was  unequiv¬ 
ocally  negative.  As  was  often  the  case,  the  presence  of  dust  in  these 
samples  greatly  reduced  the  accuracy  and  reliability  of  the  dynamic  light 
scattering  results,  especially  ac  the  lower  angles. 

Summary  of  DNA  Damage  Studies 

Tables  1,  2,  and  3  of  this  report  and  the  Summary  Table  of  our  previous 
report  (34)  show  that  our  data  overwhelmingly  indicate  no  detectable  damage 
in  irradiated  DNAs  under  any  conditions  that  we  have  studied.  The  sole  ex¬ 
periment  that  indicated  positive  breakage  has  now  been  superceded  by  numerous 
experiments  at  much  higher  energy  densities,  all  of  which  give  unequivocal 
negative  results.  Therefore,  that  single  result  is  now  presumed  to  be 
spurious. 

Numerous  direct  attempts  to  produce  a  damaged  DNA  similar  to  that  ob¬ 
served  in  the  spurious  experiment  by  conceivable  experimental  errors,  such 
as  deliberately  leaving  residual  chromic  acid  cleaning  solution  in  the  radi¬ 
ation  cell,  were  unsuccessful.  However,  in  the  course  of  subsequent  noncon¬ 
tract  research,  we  have  actually  twice  experienced  similar  damage  in  DNAs 
stored  in  dilute  solution  in  our  refrigerator,  presumably  a  consequence  of 
the  rapid  production  of  nucleases  by  contaminating  cold-insensitive  micro¬ 
organisms.  We  are,  therefore,  strongly  inclined  to  ignore  that  one  experi¬ 
ment,  and  accept  the  obvious  conclusion  that  irradiated  DNAs  suffer  no  dam¬ 
age  detectable  by  our  techniques  at  the  energy  densities  attainable. 
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TABLE  3.  SUMMARY  OF  TWO-WAVELENGTH  DNA  IRRADIATION  EXPERIMENTS 

(1060  nm  +  530  nm) 


ATTEMPTS  TO  OBSERVE  TRANSIENT  DISTORTION  BIREFRINGENCE  OF  DNA 
INDUCED  BY  PICOSECOND  PULSES  OF  1060-nm  LASER  LIGHT 


One  of  the  potential  damage  mechanisms  for  DNA  conjectured  in  the  previ¬ 
ous  report  (34)  was  based  on  intramolecular  strain  caused  by  light-induced 
torques  on  the  optically  anisotropic  DNA  bases.  Even  though  strain  to  the 
point  of  fragmentation  did  not  occur,  strain  may  still  take  place  to  a  suffi¬ 
cient  extent  to  produce  a  detectable  birefringence.  An  experiment  to  test 
this  hypothesis  is  worthwhile  for  several  reasons.  If  observed,  such  an  in¬ 
tramolecular  strain  might  permit  a  rough  estimate  of  the  laser  powers  ulti¬ 
mately  required  for  fragmentation.  Also,  if  observed,  the  transient  decay  of 
the  optically  induced  birefringence  could  provide  new  information  regarding 
the  dynamics  of  these  macromolecules  on  the  picosecond  time  scale.  Even  if 
not  observed,  it  might  be  possible  to  set  a  useful  lower  limit  on  the  re¬ 
straining  forces  governing  tilting  of  the  DNA  bases  with  respect  to  the  helix 
axis,  about  which  essentially  nothing  is  known. 

Theory  of  the  Optical  Kerr  Effect  in  Liquids  Comprised 
of  Single  Rigid  Molecules 

In  a  liquid  comprised  of  cyl indrically  symmetric  nonpolar  molecules  with 
Dolarizability  aQ  =  +  2a,  )/3  and  anisotropy  Aa  =  the  dipole 

moment  (uzi)  of  the  ith  molecule  induced  by  a  z-polarized  constant  macroscopic 
electric  field  (Ez)  can  be  written  in  the  form 

u z i  =  fEzl>oPo(cos  'V  +  I  Aa  P2^cos  e i ) ] »  (1) 

where  n..  is  the  polar  angle  between  the  molecular  symmetry  axis  and  the  electric 
field,  P  (cos  ■•.)  is  the  ?th  Legendre  polynomial  of  cos  ^ ,  and  f  =  E  -j  Qc/Ez  is 
the  internal  field  correction  factor  relating  the  local  cavity  field  E^Qc  to 


Interaction  energy  between  the  applied  electric  field  Ez°  =  rEz  and  the 
total  induced  polarization  in  the  sample  of  N  molecules  in  volume  V  is  (41) 


Neglecting  the  orientation  independent  part,  and  generalizing  to  the  mean  in¬ 
teraction  energy  with  a  very  rapidly  oscillating  optical  field  (for  which 
(E/7/2  -»  (E2°)2/4),  this  relation  takes  the  canonical  form  of  Kubo  (42): 

H(t)  =  - F ( t )  •  A  =  - 
where  the  second  equality  serves  to  define  the  quantities  F(t)  and  A.  The 

o  p 

time-dependence  of  F(t)  =  (E  )  /4  is  assumed  to  be  very  slow  compared  to  the 
optical  frequency  of  the  irradiating  light. 

2 

Using  the  relations  Dz  =  Ez  +  4tiPz  =  czE2>  nz  =  ez  (at  optical  frequen¬ 
cies),  and  nz  =  nQ  +  An  ,  it  is  readily  found  to  lowest  order  in  A a  that 


(Ezu) 


f[|  Aa  J  P2  (COS  0 j )] , 


(3) 


$n. 


4tt  f  Aa 

3~  n  V 
o 


[  p2  (cos  e.)  =  B  , 


(4) 


where  the  last  equality  serves  to  define  B  =  6nz> 

From  Kubo's  fluctuation-dissipation  theory  (42),  the  linear  response  of 
B  -  cnz  to  a  sinusoidally  varying  F(t)  =  [ (Ez°)?/4)cos  wt  (where  u>  is  here 
understood  to  be  much  smaller  than  optical  frequencies)  is  simply 


B(t)  =  Re  BA ^ } Fo  cos  at  +  Imf  ’BA^^Fo  Sin  wt  ’  (5) 

where  the  real  and  imaginary  parts  of  the  susceptibility  are  given  by 


Re{>BA(w)  ‘  -  py-  f  <AB(t)>  sin  ut  dt 

B  B  J 

o 


(6) 


Imf 


,-:T  1  --AB(t)'  cos  „t  dt. 

V  J 

r. 


(7) 


Here  T  is  the  absolute  temperature,  k^  is  Boltzmann's  constant,  and  the  angu¬ 
lar  brackets  denote  equilibrium  averages  for  the  system  in  the  absence  of  the 
perturbing  optical  field. 
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1^^ 
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For  f(t)  varying  slowly  compared  to  the  rotational  relaxation  time  ( e . g . , 
of  CS2  in  the  present  instance).  Inn  -  ^( .  ) }  =  0,  and  Re  {/g^(ui) )  -  AB>/kgT. 
Thus  the  "steady  state"  (i.e.,  . =0)  response  of  the  refractive  index  to  the 
perturbing  optical  pulse  is 


(E. 


nn  = 
z 


kBT 


8r 
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X  ( cos  ■)i)  P 2 ( cos  j)  ■. 
i  j 


(8) 


Now,  using  the  identity  P2(cos  •  i )  =  /  “ef  Y2g(si. ) ,  where  ft.  =  (0.,*..),  and 
Y  (r..)  is  the  corresponding  spherical  harmonic,  one  obtains 


P2(cos  i ) 


P2(cos 


J- 


ill  <y 

5  '20 


(fii)Y20(ftj)> 


1  <P2<Vi 


(9) 


where  the  last  line  follows  from  the  addition  theorem  of  spherical  harmonics 
(43)  and  the  fact  that  the  liquid  exhibits  an  isotropic  equilibrium  state. 

The  .  j  are  unit  vectors  directed  along  the  symmetry  axes  of  molecules  i 
and  j. 

The  instantaneous  angular  distribution  of  symmetry  axes  can  always  be 
expanded  in  the  complete  set  of  Legendre  polynomial s .  It  is  apparent  from 
Eq.  4  that  6nz  must  be  proportional  to  the  coefficient  c2(t)  of  the  P2(cos  • ) 
term,  provided  that  c 2 ( t )  vanishes  in  the  equilibrium  state  and  does  not  con¬ 
tribute  to  n  ,  as  in  the  case  for  isotropic  fluids.  Under  the  same  conditions, 
still  to  lowest  order  in  ;.a,  it  is  readily  shown  that  6n  =  :,n  =  -2*n  . 

Thus,  one  obtains  finally 

An  =  6n  -  6n  =  3An 

Z  X  L. 

(10) 

Mrs)  S2  2  .  n  2 

rTkTT  ~T5  {1  +  (N_1)  ‘  P2(l‘i  *  ^Ez  ^  ’ 

0  D 

where  (Ez°)2  is  the  square  of  the  (peak)  amplitude  of  the  perturbing  optical 
beam  ijn  vacuo.  This  formula  differs  from  the  standard  formula  applicable  to 
dilute  gas  phase  molecules, 
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15  n  kDT 
o  B 


(Ez°)2, 


2n(g)  Aa2 

15n  kRT 

0  D 


(Ez°)‘ 


(ID 


by  the  appearance  of  the  factor  in  square  brackets.  In  Eq.  11,  C  is  the 

3  9  q 

concentration  in  g/cm  ,  Ag  =  A«/V  ,  where  V  =  molecular  volume  in  cm  ,  and 

hi  m 

V  =  specific  volume  (cnr  per  g);  hence  Cg(V/Vm)  =  N/V. 

The  Clausi us-Mosotti  theory  of  the  internal  field  gives 


_3  ,V,  n02  '  1 

‘o  4W  ~~2  7 

nfl  +  2 


_3  no  ~  1 

4n  (navp  n  4  +  2 
o 


(12) 


where  M,  Nftv,  and  p  are  the  molecular  weight  in  daltons,  Avogardro’s  number, 
and  the  density,  respectively.  The  same  internal  field  theory,  which  is 
known  to  be  only  approximate,  gives 


f 


(13) 


Unfortunately,  the  value  of  the  orientation  pair  correlation  function 
0PC2  =  fl  +  (N  -  1 ) <^2( u i  ’  ]i?)> J  1S  not  available  for  CS2,  and  no  simple 
theory  is  available  for  guidance.  Values  reported  for  chloroform  and  nitro¬ 
benzene  are  2.0  and  2.8,  respectively. 


Estimation  of  Molecular  Parameters  for  CS2 

For  CS„  at  25°C,  the  Handbook  of  Chemistry  and  Physics  gives  a  refractive 
c  3 

index  n  =  1.628  and  a  density  p  =  1.595  q/cm  .  Using  these  values  in  Eq.  12 
0  -29  3 

above  gives  «  =  6.7  x  10  cm  .  The  value  f  =  1.55  is  estimated  from  Eq.  13. 

The  birefringence  induced  by  high-power  optical  pulses  has  been  found 
experimentally  to  obey  the  relation  (44) 
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-  ?0  '> 

An  =  2.2  x  10  •  (vol  ts/ia)RMS 

=  2.0  x  10  11  (statvolts/ciii)RR|<. 

=  1-0  x  10-11  (Ez0)2.  (14) 


Comparing  this  expression  with  Eq.  10  and  using  the  values  of  nQ,  a  ,  f,  and 
N/V 


iN^v/M  given  above,  yields  finally 


A  a  = 


7.25  x  10 
(0PC2) 


■24 


cm 


(15) 


The  interesting  difficulty  here  is  that  if  0PC2  were  of  order  1.0,  then 

tn  would  be  comparable  to,  or  even  larger  than,  the  estimated  a  =  6.7  x 

-24  3  0 

10  cm  .  A  second  estimate  of  a  is  available  from  the  molar  polarization 

3  ® 

(4rt/3)N»  a  =  21.1  cm  /mole  of  CS„  dissolved  in  other  nonpolar  solvents,  as 

MV  0  ^  _p4  _3 

tabulated  by  Debye  (45).  This  datum  gives  aQ  =  8.4  x  10  cm  ,  comparable 
to  the  previous  value.  One  is  forced  to  conclude  that  either  Aa  or  0PC2,  or 
both,  are  much  larger  than  expected.  Typically  Aa  =  is  only  a  few 

percent  of  a  “  (l/3)(-i^  +  2a  ) .  However,  for  CS,,,  we  must  have  Aa/«0 
greater  than  0.3  if  0PC2  is  less  than  10! 

The  point  of  the  above  discussion  is  to  note  that  CS^  provides  an  extra¬ 
ordinarily  large  induced  bi refri ngence  due  to  either  an  enormously  large  Aa 
in  comparison  to  n  or  to  a  very  large  equilibrium  orientation  pair  correla¬ 
tion  function,  which  implies  that  many  CS^  molecules  will  orient  cooperatively 
due  to  strong  intrinsic  orientational  coupling  among  them,  or  to  both  effects. 
From  a  comparison  of  Eqs.  10  and  11,  we  see  that  the  simple  standard  formula 
in  Eq.  11  can  be  used  if  CS^  is  assigned  an  apparent  Aagp  defined  by 


app 


41.55)  (0PC2) : 


,  9  v  in-23  3 

1.2  x  10  cm  . 
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(16) 


where  use  has  been  made  of  Eq.  15.  This  apparent  anisotropy  substantially 

exceeds  a  for  the  reasons  cited  above, 
o 

Theory  of  Induced  Distortion  Birefringence  of  DNA 

An  optical  pulse  of  duration  10"^  sec  passes  by  too  rapidly  to  effect 
significant  reorientation  of  the  helix-axis  of  DNA  but  may  be  able  to  tilt 
the  base-pairs  (which  are  regarded  here  as  a  single  rigid  unit)  with  respect 
to  the  helix-axis,  thereby  inducing  some  distortion  birefringence. 

Attention  is  focused  on  a  single  base-pair,  which  is  imagined  to  have  an 
optical  polarizability  that  is  cylindrically  symmetric  about  the  fixed  local 
helix-axis  in  the  undistorted  configuration.  That  fixed  helix-axis  in  turn 
makes  an  angle  e  with  respect  to  the  polarization  vector  z  of  the  optical 
pulse.  If  there  were  no  distortion,  the  interaction  energy  of  the  single 
base-pair  with  the  optical  field  would  be  given  by  the  appropriate  modifica¬ 
tion  of  Eq.  2  for  a  single  molecule. 

The  orientation  of  the  base-pair  with  respect  to  the  helix-axis  is  gov- 

2 

erned  by  an  angular  harmonic  restoring  potential  gn  /2,  where  g  is  the  restor¬ 
ing  torque  constant  and  n  is  the  angle  of  deflection  or  tilt  away  from  the 
eauilibrium  configuration.  The  equilibrium  probability  distribution  for  n  is 

2 

e  2kBT 

P(n)dn  =  ■  -  dr,.  (17) 

/2"kBT/g 


The  interaction  energy  between  the  optical  pulse  and  the  base-pair  is 


(E  °)2 

H(t)  = - f{a0Po  +  !  P2  [cos  (0  "  n)]} 


(18) 


If  g  is  sufficiently  large,  then  n  -  1.  Direct  expansion  for  small  n  1 

gives 

2  2 

P2[cos(o  -  n)]  -  (l-2r  )P2[cos  o]  +  -7J-  +  Jr  cos  u  sin  o  .  (19) 
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Because  the  helix-axes  are  held  fixed,  the  system  can  respond  to  the  optical 
pulse  only  by  changing  r,.  Thus,  the  interaction  energy  between  the  optical 
pulse  and  those  degrees  of  freedom  (i.e.,  n)  that  are  free  to  respond  is 


H  (t)  =  - 


f(Ez  l  !  |2a,x  2  2  ^ 

4  |  j  3  '  i_2'i  r^Lcos  u]  +  n  +  3r.  cos  sin  ;/ j 

F(t) -An  . 


(20) 


The  change  in  refractive  index  induced  by  the  pulse  acting  on  this  one 
base- pair  is  found  from  Eq.  4  in  the  form 


6nz  =  T  n  r  P2  [cos  (0  '  (2D 

o 

Because  the  a  value  for  a  given  base-pair  is  fixed,  only  that  part  of  <5n 
that  varies  with  n  will  actually  contribute  to  a  6n  induced  by  base-tilting. 
Therefore  we  keep  only  the  lowest  n-dependent  terms 


4  ;-  f  2  r  _ 

Sn2  ~  T~  t-2ri  P2^cos  °)  +  +  3n  sin  e  cos  6} 


=  B  . 

n 


(22) 


Again,  linear  response  theory  may  be  employed  with  the  result  that  the 
steady-state  birefringence  is  given  by 


r(Ez°)21 


<A  B  : 

n  n 

kBT 


(23) 


(Ez  ^  f 2  8tt  a*2  .  -  2  n2  2 

4~  9 - r  <t_2n  p2‘-cos  +  +  3n  sin  0  cos  }  >5 


where  the  angular  brackets  denote  an  average  over  n  using  the  canonical  dis¬ 
tribution  in  Eq.  17.  Odd  powers  of  n  average  to  zero,  and  fourth-order  terms 

2 

are  neglected;  so  only  the  -'n  =  k^T/g  •  1  terms  contribute.  It  is  assumed 

that  we  are  in  the  stiff  limit  where  g  :■>  kDT.  Consideration  of  a  collection 

D 


27 


of  N  independently  tilting  base-pairs  can  be  accommodated  by  multiplying  Eq. 
23  by  N  and  performing  an  angular  average  over  the  uniform  distribution  of 

helix-c  is  orientations  e.  The  entire  final  contribution  comes  from  the 

2 

< ( 3 n  sin  cos  i)  >  term,  and  is  given  by 


c_  _  4i;  fc  Au‘  ,N,  ,r  0,2 

6nz  "  15  1  ~7T  (V}  (Ez  }  ‘ 
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Again,  An  -  -n,  -  an  =  3..n7,  so  finally 

2  2  2 

,n  _  4::  f  Act  ,N,  0. 

An  -  r  1  TT  ¥  (Ez  }  ’ 


(24) 


(25) 


which  may  be  compared  to  Eqs.  10  and  11  for  freely  orienting  species.  It 
should  be  noted  that  Eqs.  24  and  25  apply  only  when  g/k^T  »  1.0. 


Attempts  to  Observe  Transient  Distortion  Birefringence  of  DNA 

Herring  sperm  DNA  (Na  salt).  Type  VII  from  Sigma  Lot  #670-00251,  was 

sheared  and  concentrated  down  to  approximately  63  mg/ml  using  a  french  press 

at  12,000  psi,  and  then  evaporating  in  a  vacuum.  The  size  of  the  sheared 

fragments  was  estimated  to  be  roughly  200,000  daltons  with  an  average  diffu- 

Q  o 

sion  coefficient  D  *  8.2  x  10  cm  /s.  Shearing  was  necessary  to  achieve  such 
a  highly  concentrated  solution.  This  was  then  used  as  the  active  medium  in  an 
optical  Kerr  effect  streak  shutter  (46),  as  indicated  in  Figure  1.  Unfortu¬ 
nately,  no  detectable  birefringence  of  the  DNA  was  ever  observed  with  DNA  as 
the  active  medium  in  this  apparatus. 

The  intensity  of  the  maximum  transmitted  signal  appeared  as  a  vertical 
deflection  of  5  cm  on  the  oscilloscope  display.  A  deflection  as  small  as  2  mm 
would  have  been  readily  observed.  In  addition,  the  probe  beam  intensity  was 
increased  60-fold  for  the  DNA  experiment.  Thus,  the  transmission  sensitivity 
of  the  streak  shutter  was  GOx50/2  =  1500  times  higher  in  the  DNA  experiment 

than  for  CS9.  Because  the  particular  apparatus  used  in  this  streak  shutter 

c  2 
gives  a  transmitted  intensity  signal  proportional  to  An  ,  we  may  conclude 

that 
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(26) 
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n  o 

From  the  measured  intensity  of  the  gating  pulse  ( I m  10  W/cm),  the  RMS  volts/m 

can  be  recalculated  and  used  in  Eq.  14  to  give  An™?x  =  8x10*^.  Thus, 

. „ ma  x  , 

An0NA  -  2xl0'6. 


A  lower  limit  on  the  torque  constant  g  can  be  obtained  from  Eq.  26  by 
using  Eqs.  11  and  16  for  An^x  on  the  right-hand  side,  and  Eq.  24  for  An™j^  on 
the  left.  That  is. 


4-r  f2  /L2  /pNAVwP  o,2  _ 1_  2tt 

5  9  "0  m  2  -  38-7  15  kBT 


A  2  Awr  0,2 
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(27) 


where  the  tilde  (~)  denotes  quantities  pertinent  to  the  DNA  solution.  The 

narameters  for  CS2  were  given  earlier.  For  water,  no  =  1.334,  and  Eq.  13 

gives  f  =  1.26.  For  the  DNA  solution,  r  =  0.063  p/1 ,  and  M  =  662  per  base- 

_2fl~  3 

oair.  For  one  base-pair  we  take  Aa  =  -14.1x10  'erf  ,  which  we  have  calculated 
from  the  flow  birefringence  data  of  Harrington  (47)  in  0.2  M  NaCl ,  using  his 

o 

Eq.  2  together  with  the  persistence  a  =  490  A  reported  by  Voorduow  et  al .  (48). 
This  value  is  slightly  higher  than  that  (Aa  =  -12.9xlC_  enf- )  computed  by 

o 

Harrington  using  a  =  660  A.  Using  these  values,  we  obtain ,  final ly,  the  lower 
1  ini  t 


g  -  (3.23)  kgT  =  1.32xl0~13  dyne-cm.  (28) 

This  lower  limit  only  marginally  satisfies  the  inequality  g  >>  kgT  required 
for  validity  of  Eqs.  20-25,  and  in  that  sense  is  disappointingly  small.  The 
torsion  constant  between  base-pairs  for  the  twisting  of  DNA  has  recently  been 
measured  (49)  and  found  to  be  3.8x10  dyne-cm.  Thus,  this  experiment  has 
established  only  that  the  torque  constant  for  base-tilting  must  exceed  approx¬ 
imately  1/30  of  the  torsion  constant  for  twisting  successive  base-pairs 
(about  the  the  helix  axis)  with  respect  to  one  another.  That  is  neither  a 
surprising  nor  a  very  useful  result. 

The  same  DNA  preparation  likewise  showed  no  induced  birefringence  after 
being  immersed  in  boiling  water  for  20  min  (to  denature  the  DNA)  and  then 
quenched  in  an  ice  bath. 
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We  note  that  if  the  intensity  of  the  orientin']  optic.il  pulse  could  be 
increased  by  a  factor  of  10,  then  even  a  negative  result  would  provide  an 
interesting  lower  limit  for  the  base-tilting  torque  constant. 

DAMAGE  STUDIES  OF  LIPID  BILAVER  VESICLES 

Dipalmitoyl  Phosphatidyl  Choline  (DPPC) 

This  synthetic  pure  lipid  (i.e.,  DPPC)  forms  bilayer  vesicle  membrane 

structures  by  appropriate  preparation,  or  introduction,  in  aqueous  solution. 

It  does  exhibit  a  thermal  order-disorder  transition  in  the  bilayer  near 

T  -  41 °C.  The  experiments  reported  here  refer  to  vesicles  at  T  <  T  ,  which 
in  m 

exhibit  the  more  ordered,  less  fluid  phase  of  the  membrane  bilayer.  They  are, 
therefore,  probably  less  representati ve  of  biological  membranes,  which  are 
composed  of  a  variety  of  mixea  lipids  for  which  the  fluid  phase  is  believed 
to  prevail  at  room  temperature.  A  few  experiments  performed  subsequent  to  the 
contract  indicate  that  the  stabi  ity  of  DPPC  vesicle  suspensions  is  consider¬ 
ably  reduced  for  T  >  T  . 


Preparation  of  DPPC  Vesicles 

Vesicles  were  prepared  from  a  95?  ethanol  solution  of  dipalmitoyl  phos¬ 
phatidyl  choline  (DPPC)  at  a  concentration  of  23  umole/ml  (50).  A  0.50-ml 
aliquot  of  the  ethanol ic  lipid  solution  was  injected  into  10  ml  of  buffer 
(0.1  m  NaCl  and  0.01  M  Tris-HCl)  at  a  rate  of  0.20  ml/min  to  give  a  final 
lipid  concentration  of  2.3  nmol  e/ml.  This  injection  method  produced  vesicles 
about  400  A  in  radius,  in  good  agreement  with  literature  values  (51).  Some 
polydi spersi ty  existed  even  before  irradiation,  as  evidenced  by  a  decline  in 

the  apparent  diffusion  constant  at  low  scattering  angles  (e  <  45‘  ,  or 
2  10  -2 

K  <  1.0  x  10  cm  ) ,  as  shown  in  Figure  2.  After  exposure  to  picosecond 
laser  radiation,  the  2.3  umole/ml  sample  solution  was  transferred  directly  in¬ 
to  a  light  scattering  cell  that  had  been  previously  rinsed  free  of  dust  by 
washing  with  filtered  (0.2  um  Mi  1 1 i pore)  buffer.  The  sample  solution  was 
diluted  with  filtered  (0.2  um)  buffer  to  0.23  umole/ml  for  use  in  dynamic 
light  scattering  (DLS). 
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D  vs.  K  plot  for  a  sample  with  typical  pulse  train  irradiation  of  600  mJ/cm  total 
energy  density.  Laser  wavelength  for  dynamic  light  scattering  is  632.8  nm. 

(Error  bars  indicate  full  range  of  measured  D  values.) 


Experimental  Measurements  on  DPPC  Vesicles 


Nine  successful  exposures  of  the  vesicles  were  carried  out  ranginq  from 

2  2 

an  energy  density  of  9  mj/cm  in  a  single  pulse  to  3  J/cm  in  an  entire  pulse 
train.  The  exposure  parameters  of  the  various  samples  are  given  in  Table  4. 

TABLE  4.  SUMMARY  OF  DPPC  VESICLE  IRRADIATION  EXPERIMENTS 


Sample 

Irradiation 

Energy 

density 

2 

(mJ/cm  ) 

Damage 

1 

pulse  train 

706 

yes 

2 

1  pulse 

9.0 

(massive) 

possibly 

3 

pulse  train 

580 

(si ight) 
yes 

4 

pulse  train 

614 

no 

5 

pulse  train 

3,080 

yes 

6 

pulse  train 

529 

(slight) 

no 

7 

pulse  train 

533 

ambiguous 

8 

pulse  train 

533 

(probably 

no 

9 

pulse  train 

529 

no 

Studies  with  DLS  were  performed  on  all  of  the  samples  as  well  as  their 
appropriate  controls.  For  a  polydisperse  solution  of  scatterers  with  dimen¬ 
sions  small  compared  to  K”^,  where  K  =  (4m/x)sin  e/2  is  the  scattering  vec¬ 
tor,  n  the  index  of  refraction,  and  e  the  scattering  angle.  The  apparent  dif 
fusion  coefficient  obtained  from  the  dynamic  light  scattering  is  the  so- 
called  z-average  value  (52)  defined  by 


2  2 
>  n.m.  D./)  n .m . 

|  i  i  1  j  J  J 


(29) 


where  n^ ,  rru ,  and  D^  are  the  number  concentration ,  mass,  and  diffusion  coeffi 
cient  of  the  ith  species.  For  spherical  particles  of  dimension  comparable  to 
or  larger  than  K-*,  the  appropriate  general ization  of  Eg.  29  is 


(30) 
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where  P.(K)  is  the  internal  interference  factor  characterizing  the  ith  species. 
P. (K)  decreases  from  a  maximum  value  of  1.0  at  small  K  to  a  lower  value  (de¬ 
pending  on  the  particular  shape  and  size  of  the  scatterer)  at  any  given  larger 
2 

value  of  K  .  For  sufficiently  large  particles,  P.(K)  is  negligible  at  larger 

o  n  1 

K*-,  but  at  small  K  may  permit  the  large  particles  to  dominate  the  scattering 

from  the  smaller  species,  thus  causing  a  decrease  in  the  apparent  diffusion 

2 

coefficient  at  small  K  (as  indicated  in  Fig.  2).  Conversely,  a  significant 

2 

decline  in  D  at  small  K  implies  the  existence  of  very  large  scatterers  in  the 
preparation. 

The  apparent  Stokes  radius  R^  =  kgT/67rnD--where  kg  is  Boltzmann's  constant, 

T  the  absolute  temperature,  n  the  solution  viscosity,  and  D  the  diffusion  coef- 

ficient--was  computed  for  the  very  large  particles  from  the  low-angle  (25°  or 

30°)  data  for  comparison  with  static  light  scattering  measurements  of  the  radii 

of  gyration  R  of  those  same  particles.  The  latter  values  were  computed  using 
9  2  2 

the  relation  applicable  for  k  Rg  $  1.0 

iN(e)/iN;o)  =  1  -  K2Rg2/3,  (31) 

where  IN(o)  is  the  "normalized"  intensity  (i.e.,  corrected  for  the  relative 
size  of  scattering  volume  subtended  at  o),  and  1^(0)  is  the  extrapolated 
value  of  IM(o)  at  o  =  0.  A  comparison  of  Rh  and  values  obtained  for  either 
radiated  samples  or  controls  over  the  mid-  to  low-angle  range  of  measurements 
shows  that  Rg  >  Rh  (typically,  Rg  =  1,600  A,  Rh  =  720  ft) ,  which  implies  a 
nonspherical ,  irregular  or  anisometric  shape  for  particles  of  very  large  di¬ 
mension,  thus  indicating  that  they  are  not  simply  very  large  spherical 
vesicles . 

O 

The  ordinary  vesicles  (R.  z  450  A)  do  not  satisfy  the  criterion, 

2  2  ^ 

K  R^/3  >  1,  and  hence  do  not  produce  sufficient  variation  in  I^(  )  with  to 

yield  a  meaningful  estimation  of  Rg  using  the  present  data.  In  experiments 

subsequent  to  the  contract  period,  we  have  been  able  to  obtain  very  precise 

2  2  10 

static  and  dynamic  light  scattering  data  over  the  high  K  range  (K  1x10 
cm  2 ;  45°  <  )  <  120°)  for  somewhat  larger  vesicles,  and  have  established  that 
both  and  R  are  consistent  (within  a  few  ft)  with  a  spherical  shell  of  outer 
diameter  580  A  and  thickness  37  ft,  the  latter  of  which  is  the  generally  accept¬ 
ed  value  for  lipid  bilayers.  The  point  of  the  latter  remark  is  that  we  have 
now  some  hard  evidence,  for  at  least  one  preparation,  that  the  actual  shape 
of  the  ordinary  vesicle  species  is  close  to  spherical,  as  is  widely  believed. 
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Slow  secular  changes  in  the  properties  of  all  samples  (including  the 
nonirradiated  controls)  with  time,  required  that  the  samples  and  controls  be 
exposed  and  analyzed  the  very  same  day. 

Resul ts 

The  criterion  for  assessing  damage  to  the  vesicles  was  whether  the 
2 

D  vs  K  plots  for  the  irradiated  samples  changed  significantly  in  any  respect 
when  compared  with  their  respective  controls,  which  had  not  been  subjected  to 
picosecond  laser  radiation  but  otherwise  had  received  identical  handling,  in¬ 
cluding  loading  into  the  irradiation  cell.  Both  were  observed  at  the  same 
time. 

2 

Two  regions  of  the  D  vs  K  curve  are  of  special  interest.  The  high-angle 
2 

(n  •  45  ) ,  or  large  K  ,  D-values  give  an  indication  of  the  z-average  inverse 

2 

hydrodynamic  radii  of  the  ordinary  vesicles.  The  low-angle,  or  small  K  , 

D-values  contain  primarily  information  about  scattering  centers  of  larger 

dimension,  such  as  clumps  of  vesicles  or  large  lamellar  aggregates  (liposomes) 

that  scatter  light  predominantly  in  the  forward  direction. 

Exposure  of  the  vesicle  suspensions  to  entire  pulse  trains  with  energy 

densities  of  600  mJ/cm  (peak  electric  field  "..5.8x10  V/cm  )  or  greater  did 

2 

produce  a  distinct  change  in  the  D  vs  K  curve,  as  indicated  in  Figure  2-- 
which  shows  that  the  apparent  diffusion  coefficient  decreased  by  more  than  25"' 
at  ■  =25  and  by  a  smaller  but  still  significant  amount  at  higher  angles. 

Our  interpretation  of  this  result  is  that  the  number  and/or  dimensions  of  the 
nonspherical  aggregates,  presumably  multilamellar  liposomes,  have  increased  at 
the  expense  of  the  smaller  ordinary  vesicles  destroyed  by  the  light  pulse. 

Data  for  the  various  samples  are  summarized  in  a  general  way  in  Table  4. 
Values  of  the  apparent  diffusion  coefficient  observed  at  various  angles  are 
recorded  in  Table  5. 

The  massive  damage  sustained  in  sample  1,  which  was  exposed  to  an  entire 
pulse  train,  completely  overshadowed  the  more  modest  changes  in  sample  2, 
which  was  exposed  t ^  a  single  pulse.  However,  after  observing  the  behavior 
of  the  remaining  samples,  which  include  some  samples  of  clear-cut  negative 
results  as  well  as  some  significant  examples  of  more  modest  damage,  including 
sample  3  (shown  in  Fig.  2),  it  seems  possible  that  sample  2  actually  sustained 
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TABLE  5.  APPARENT  DIFFUSION  COEFFICIENTS  (DxlO°  cnT/sec)  OF  DPPC  VESICLES 
DETERMINED  BY  DYNAMIC  LIGHT  SCATTERING 


Control  for  8  &  9  3.44  3.97  4.70  5.09  5.43  5.67 


2 

some  significant  damage  after  all.  The  pulse  energy  of  9  mJ/om  in  that 
single  pulse  corresponded  to  a  peak  electric-field  strength  (assumed  uniform 
over  a  10-psec  pulse)  of  ^7. 1x10  V/cm. 

Curiously,  sample  5,  which  was  exposed  to  an  entire  pulse  train  at  very 

high  energy  density,  showed  rather  slight  indication  of  damage  in  the  sense 

?  ? 

that  there  was  a  decrease  in  D  only  at  small  K  ,  none  at  large  K  .  However, 
that  difference  between  the  irradiated  sample  and  control  persisted  into  the 
second  day  after  exposure,  as  indicated  in  Table  6,  thereby  increasing  our 
confidence  in  that  result.  That  sample  was  irradiated  in  the  longer,  narrower 
bore  cell  to  achieve  the  higher  energy  density.  These  vesicles  are  known  to 
stick  to  glass,  and  the  large  increase  in  surface  area  may  have  "fractionated" 
the  irradiated  sample  in  an  unknown  way. 

TABLE  6.  APPARENT  DIFFUSION  COEFFICIENTS  OF  DPPC  VESICLES 
AS  FUNCTIONS  OF  SAMPLE  AGE 


(0  x  10^  cm^/sec) 

Sample 

Age 

6  =  25°  9  =  90° 

CD 

It 

1 — ‘ 
r\i 
o 

° 

Damage 

3 

day  1 

4.53 

7.29 

yes 

day  2 

1.49 

6.90 

4 

day  1 

5.82 

7.55 

no 

day  2 

4.97 

7.22 

5 

day  1 

1.40  5.07 

si ight 

day  2 

0.932  4.72 

Control 

day  1 

2.09  5.01 

for  5 

day  2 

1.21  5.12 

2 

These  results  indicate  an  approximate  damage  threshold  of  about  600  mj/cm 
for  DPPC  vesicles  exposed  to  trains  of  picosecond  pulses  of  1060-nm  radiation. 
In  addition,  there  is  some  indication  that  more  modest  damage  might  be  inflict- 
ed  by  single  picosecond  pulses  with  an  energy  density  of  9  mj/cm  .  It  is 
possible  that  electrostrictive  or  other  forces  present  during  the  pulse  actual¬ 
ly  rupture  the  vesicle  every  time,  but  that  substantial  recovery  occurs  be¬ 
tween  pulses  in  these  systems.  Certainly,  vesicle  fragments  do  not  move  far 
-9 

in  10  sec,  and  lipids  are  insoluble  in  aqueous  solution.  Thus,  the  observed 
damage  may  be  minimal  unless  the  vesicle  is  subjected  to  repeated  battering 
by  closely  spaced  pulses,  as  in  a  pulse  train.  In  a  system  where  chemical 
potential  gradients  and  osmotic  pressure  differences  exist  across  membrane,  as 
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in  living  systems,  a  single  rupture  might  lead  to  complete  destruction  of  the 
osmotically  strained  vesicle,  so  the  damage  threshold  in  terms  of  peak  elec¬ 
tric  field  might  be  the  same  for  single  pulses  and  whole  pulse  trains  in  that 
case.  Indeed,  our  studies  of  damage  thresholds  in  primate  retinae  support 
this  hypothesis.  (See  Part  II  of  this  report.) 

It  is  also  evident  from  Table  6  that  the  differences  between  samples  3 
(damaged)  and  4  (undamaged)  persist  into  the  second  day  after  exposure,  despite 
the  long-term  changes  in  both  samples.  Again,  this  argues  that  the  differences 
between  these  samples  are  genuine. 

The  possibility  of  a  thermal  damage  mechanism  playing  a  role  in  the  above 
experiments  can  be  unequivocally  ruled  out--the  temperature  rise  in  our  samples, 
as  estimated  from  the  difference  in  energy  between  the  incident  and  the  trans¬ 
mitted  pulse  trains,  was  less  than  0.02'C. 

Egg-Yolk  Lecithin  (EYL)  Vesicles 

Our  results  for  DPPC  vesicles  apply  to  the  ordered  low-T  phase  of  the 
membranes  prevailing  below  T^  *  41'C.  The  question  naturally  arises  whether 
similar  results  would  be  found  for  membranes  comprised  of  natural  lipids,  such 
as  egg-yolk  lecithin  (EYL).  Although  both  lipids  are  phosphatidyl  cholines, 
the  DPPC  has  two  identical  saturated  18-carbon  fatty  acid  chains  esterified  to 
the  glyceryl  moiety,  whereas  EYL  phosphatidyl  cholines  generally  contain  a 
saturated  fatty  acid  esterified  at  the  1  position  in  addition  to  the  unsatu¬ 
rated  fatty  acid  (usually  oleic)  at  the  2  position.  As  a  result  of  this  chem¬ 
ical  difference,  and  also  the  mixture  of  phosphatidyl  cholines  present,  the 
EYL  membranes  exhibit  a  broad  thermal  transition  between  -15°  and  -7°C  and  are 
believed  to  exist  in  the  fluid  phase  at  room  temperature,  r'20°C. 

Preparation  of  EYL  Vesicles 

Considerable  experimentation  with  parameters  governing  the  vesicle 
preparations  was  conducted.  The  effects  of  mixing-temperature,  initial  con¬ 
centration,  needle  gauge,  and  injection  speed  were  all  investigated.  The 

2 

D  vs  K  curves  for  some  of  the  resulting  suspensions  are  indicated  in  Figure  3. 
The  final  procedure  is  described  below. 
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Figure  3.  D  vs.  K  for  different  preparations  of  EYL  vesicles. 

Prep.  I.  Step  1  (see  text  for  steps)  of  the  procedure  was  modified 
as  follows:  Inject  0.94  ml  of  10  mg/ml  stock  into  warm 
(40°C)  buffer.  Step  3  is  omitted. 

Prep.  II.  Step  2  was  modified  by  injecting  into  warm  (40°C)  buffer. 
Step  3  was  omitted. 

Prep.  III.  Step  2  was  modified  by  injecting  into  warm  (40°C)  buffer. 
Prep.  IV.  Step  2  was  modified  by  using  0.94  ml  of  10  mg/ml  stock. 
Prep.  V.  Same  as  step  4. 


For  EYL  commercially  supplied  in  ethanolic  solution,  the  following  steps 
are  carried  out: 

(1)  Rinse  out  a  25-ml  flask  and  flush  with  filtered  (0.22-um  Millipore) 
buffer,  leaving  the  magnetic  stirring  rod  in  the  buffer  solution  to 
be  rinsed  at  the  same  time. 

(2)  Withdraw  0.5  of  10  mg/ml  stock  solution  (95"  ethanol),  allow  it  to 
warm  up  for  several  minutes  in  syringe,  then  slowly  inject  into  a 
25-ml  flask  containing  10  ml  of  filtered  (0.22-um  Millipore)  buffer. 
Use  a  Hamilton  1-cc  syringe  with  a  22-gauge  needle.  The  concentra¬ 
tion  in  the  flask  is  now  0.47  mg/ml. 

(3)  Filter  the  suspension  through  3.0-um  Millipore  filters  twice. 

(4)  Clean  out  scattering  vessel  with  0.22-um  filtered  buffer.  Transfer 
0.45  ml  of  solution  into  cell  and  dilute  with  9  ml  of  clean  buffer. 
Final  concentration  is  0.22  mg/ml  lipid. 

For  lyophilized  EYL  make  a  10  mg/ml  solution  in  95/  ethanol  first,  then 

follow  steps  1-4  above. 

The  buffer  used  for  all  experiments  was  0.1  M  NaCl  +  0.01  M  Tris  (Sigma 

Trizma-base  MW  121.1)  titrated  with  HC1  to  pH  7.0. 

2  2 
The  D  vs  K  curves  in  Figure  3  show  no  limiting  plateau  at  high  K  ,  which 

indicates  a  greater  degree  of  polydispersity  than  found  for  DPPC  vesicles. 

Damage  Studies  of  EYL  Vesicles 

EYL  vesicles  were  exposed  to  mode-locked  1060-nm  laser  pulse  trains  with 

P 

energy  densities  between  400  and  660  mJ/cm  .  The  results  of  dynamic  light 
scattering  measurements  on  the  irradiated  samples  and  corresponding  identical¬ 
ly  treated  controls  are  presented  in  Table  7.  Although  samples  6,  7,  and  8, 

2 

irradiated  at  the  higher  energy  densities  (657,  652,  and  614  mJ/cm  respective¬ 
ly),  exhibit  D  values  marginally  smaller  than  those  of  the  controls  at  every 
angle,  the  data  are  not  sufficiently  different  to  warrant  a  definitive  conclu¬ 
sion  that  damage  occurred. 

2 

Unfortunately,  all  exposures  at  energy  densities  exceeding  600  mJ/cm 
were  obtained  by  decreasing  the  beam  diameter  by  20%,  so  not  all  of  the  sample 
solution  was  exposed  to  the  full  power  of  the  irradiating  beam.  This  practice 
was  necessitated  by  the  comparatively  low  energy  output  of  the  laser  even  with 
all  operating  parameters  optimized. 

Although  the  D  vs  K  curves  differed  only  marginally  from  the  controls 
for  the  most  intensely  exposed  samples,  it  is  still  possible  that  greater 
differences  were  manifested  in  the  time-course  of  the  decay  of  the  correlation 


TABLE  7.  VALUES  OF  D  AT  VARIOUS  K2  FOR  EYL  VESICLES 
DxlO8  (cm/sec) 


Sample 

E 

2 

(mJ/cm 

K2  = 

0,4 7x10 10 

K2  = 

1 . 03xl010 

k2  = 

1 , 75xlQ10 

K2  = 

3.51xl010 

k2  = 

5 . 26xl010 

1 

400 

3.68 

4.32 

4.89 

5.72 

6.14 

2 

420 

3.51 

4.03 

5.07 

5.79 

6.18 

Control 

3.92 

4.40 

5.11 

5.81 

6.22 

3 

420 

3.35 

4.40 

4.97 

5.45 

5.91 

4 

516 

3.05 

4.17 

4.94 

5.48 

5.88 

Control 

3.51 

4.16 

4.73 

5.60 

6.03 

5 

434 

3.20 

3.89 

4.49 

5.03 

5.48 

6 

651 

3.10 

3.82 

4.35 

4.88 

5.29 

Control 

3.18 

3.98 

4.39 

5.09 

5.46 

7 

652 

3.15 

4.01 

4.52 

5.01 

5.43 

8 

614 

3.02 

3.80 

4.37 

4.90 

5.27 

Control 

3.27 

4.09 

4.63 

5.07 

5.57 

function. 

That  is, 

,  the  quality 

of  single-exponential 

fit  may  still  have 

varied  significantly  between  irradiated  samples  and  controls.  A  method  for 
quantitatively  measuring  this  deviation  from  single-exponential  decay  is 
that  of  cumulant  analysis. 

The  normalized  intensity  autocorrelation  function  gv  '(t)  is  a  4th-order 

correlation  in  the  scattered  electric  field.  Under  the  assumption  of  a 

1 2) 

Gaussianly  distributed  scattered  field,  gv  '(-i)  may  be  related  to  the  2nd- 
order  correlation  function  g'^'(-r)  by  (52) 

g(2)(x)  =  1  +  6|g(1) (t) i 2  (32) 

where  p.  is  of  order  1.  If,  furthermore,  the  sample  is  polydisperse,  g^^(t) 

is  weighted  by  the  relative  intensity  distribution  of  scatterers  G(r)  with 

2 

different  decay  constants  r  =  K  D, 


g(1)(x) 


=  |  G(r)e  rxdr  e  <e 
0 


(33) 


The  exponential  may  be  written  in  the  form 


where 


r  :  |  G(r )r  dr  .  (34) 

o 

The  factor  e  1  may  be  expanded  and  inserted  in  Eq.  33  with  the  following 
result: 


9{1)(t)! 


(35) 


where 

u2  -  G(r)(r-r)2dr 

. 

u3  h  G(r)(r-r)3dr  . 

Thus,  the  intensity  autocorrelation  function  may  be  expressed  in  terms  of  the 
moments  of  G ( I ) .  It  is  then  easily  seen  that  (52) 

2  3 

^■ln[g^(x)  -  1]  =  jlrm  -  J't  +  |r  -  jp  u3  +  —  (36) 

where  use  has  been  made  of  the  expansion 

ln( 1+z)  =  z-z2/2  + 

Equation  36  is  a  polynomial  in  time  whose  coefficients  are  the  moments  of 
the  distribution  G(r).  The  polydispersity  of  a  sample  of  scatterers  is  here 
defined  as  the  ratio  of  the  second  moment  to  the  square  of  the  average,  that 
is, 

_2 

polydispersity  i^/i 

A  program  was  written  to  analyze  certain  sample  runs  according  to  the 
above  polynomial.  The  algorithm  first  subtracts  a  baseline  from  the  raw  auto¬ 
correlation  data.  This  baseline  is  first  taken  to  be  either  the  average  of 


the  last  25  points  or  the  baseline  resulting  from  a  single-exponential  fit. 

Then  the  algorithm  takes  the  natural  logarithm,  does  a  1  east-squares  fit  to 

2 

the  polynomial,  and  calculates  ■<  ,  the  sum  of  the  residuals,  which  is  compared 

2  2 
with  the  previous  x  •  The  baseline  is  varied  until  a  minimum  in  >  is  obtain- 

-  -  -1  -2 
ed.  The  corresponding  values  of  2r,  (2r)  ,  and  polydispersi ty  (u^/f  ) 

are  printed  on  the  teletype. 

The  polydispersities  obtained  for  correlation  functions  at  small  scatter¬ 
ing  angles  (9=30)  always  give  higher  values  than  at  9=120°,  as  shown  in 
Table  8.  The  high-angle  polydi spersities  were  generally  quite  similar  (i.e., 
within  10')  for  the  irradiated  samples  and  their  respective  controls,  whereas 
the  polydispersities  at  o=30"  usually  were  substantially  larger  (by  a  factor 
1.3)  for  the  irradiated  samples  than  for  controls,  as  also  shown  in  Table  8. 
Thus,  the  cumulant  analysis  indicates  an  appreciable  change  in  the  irradiated 
samples  studied  (at  435  mJ/cm  and  650  mJ/cm  )  with  respect  to  their  controls. 
Unfortunately,  this  change  was  not  consistently  manifested.  Moreover,  cumu¬ 
lant  analysis  has  some  inherent  weaknesses.  The  final  results  are  inordinately 

2 

sensitive  to  the  value  of  the  subtracted  baseline,  because  the  minimum  in  y 
with  respect  to  baseline  is  not  always  unique.  Moreover,  as  the  value  of  the 
baseline  was  varied,  the  number  of  data  points  to  be  fitted  also  changed  to 
prevent  a  zero  or  negative  argument  of  the  logarithm.  This  latter  feature  had 
the  consequence  that  the  fitting  algorithm  in  some  cases  simply  slid  progres¬ 
sively  to  higher  baselines  with  ever-fewer  data  points  without  converging  to 

2 

any  local  minimum  in  v  . 


TABLE  8.  POLYDISPERSITY  u?/f2  at  0  =  30°  AND  0  =  120 ' 

rnn  rvi  M/rcm  rc 


FOR  EYL  VESICLES 

Sample  5 

Sample  6 

435  mJ/cm2 

650  mJ/cm2 

Control 

o  =  30" 

K2  =  0.47xl010 

0.233 

0.200 

0.127 

a  =  120° 

1C  =  5 . 3xl0lu 

0.125 

0.117 

0.130 

Ratio 

(30/120) 

1  .86 

1.71 

1.02 
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Sumnary  of  EYL  Results 


The  EYL  lecithin  vesicles  proved  considerably  more  intractable  than  DPPC 
vesicles  in  every  regard  from  preparation  to  analysis  of  their  dynamic  light 
scattering  autocorrelation  functions.  Although  some  evidence  of  optical 
stress-induced  changes,  especially  in  the  polydispersity  parameter,  was  ob¬ 
served  in  some  preparations ,  we  feel  that  the  evidence  overall  is  too  slim 
to  warrant  at  this  time  an  unequivocal  statement  concerning  damage  of  these 
vesicles. 


Future  Vesicle  Work 

The  next  logical  step  is  to  examine  the  susceptibility  to  picosecond 
laser  pulses  of  osmotically  strained  vesicles;  for  example,  containing  concen¬ 
trations  of  salt  inside  differing  from  those  prevailing  outside.  In  that 
case,  which  corresponds  more  closely  to  living  organisms,  one  has  the  possi¬ 
bility  that  a  single  rupture  event  will  be  amplified  to  complete  vesicle 
destruction  by  the  osmotic  gradient,  with  less  likelihood  of  recovery.  In 
such  a  case  the  damage  threshold  for  single  pulses  could  well  approximate 
1/N  of  that  for  entire  pulse  trains,  where  N  is  the  total  number  of  pulses 
in  the  train.  In  other  words,  the  first  pulse  to  reach  the  threshold  field 
would  destroy  the  vesicles,  and  all  subsequent  pulses  in  the  train  would 
simply  add  "insult  to  injury." 
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PART  II:  PICOSECOND  OCULAR  DAMAGE  STUDIES  ON  PRIMATES 


INTRODUCTION 

This  phase  of  the  program  was  designed  to  determine  ultrashort  pulse- 
induced  ocular  damage  thresholds  in  the  primate  Macaca  fasci cul ari s*  at  two 
different  wavelengths,  530  and  265  nm,  obtained  by  frequency  upconversion  of 
the  Nd:Glass  laser  output,  and  to  compare  the  effects  of  single  pulses  and 
entire  mode-locked  pulse  trains  at  the  530-nm  wavelength.  In  the  visible 
spectrum  the  cornea,  lens,  and  vitreous  are  highly  transparent  and  thus 
damage  is  sustained  in  the  chorioretinal  region,  whereas  at  the  UV  wavelength 
the  incident  radiation  is  absorbed  mostly  within  the  corneal  epithelium. 

The  irradiation  apparatus  used  in  both  the  retinal  and  corneal  damage 
studies  was  similar  to  that  used  in  the  macromolecul ar  exposure  experiments. 
(See  Appendix  A.)  The  output  of  the  Nd:Glass  laser  was  frequency-doubled  to 
generate  530-nm  light  for  the  retinal  damage  studies  and  frequency-quadrupled 
to  generate  265-nm  radiation  for  the  corneal  threshold  experiments.  Addition¬ 
al  details  are  given  in  the  following  section,. 

The  experimental  protocol,  parameters,  and  criteria  we  used  followed  as 
closely  as  possible  those  used  in  the  bulk  of  previously  published  work,  in 
order  to  permit  direct  comparisons  with  data  acquired  in  prior  research  by 
other  investigators. 

We  describe  first  the  retinal  damage  threshold  studies,  which  comprise 
two  separate  bodies  of  experimental  data:  (1)  damage  thresholds  for  entire 
mode-locked  pulse  trains,  and  (2)  damage  thresholds  for  single  ultrashort 
pulses.  The  corneal  ultraviolet  irradiation  experiments  are  discussed  in 
the  final  section  of  this  report. 


*0riginally  we  planned  to  use  rhesus  monkeys  (Macaca  mul atta) ;  however, 
the  continuing  embargo  on  the  exportation  of  this  species  by  India  and  the 
resulting  severe  shortage  dictated  the  use  of  M.  fascicularis.  The  structure 
and  pigmentation  of  the  retinae  of  the  two  species  are  very  similar,  allowing 
direct  comparison  of  the  present  work  with  past  results. 
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RETINAL  DAMAGE  THRESHOLDS  INDUCED  BY  PICOSECOND  530-nm  LIGHT  PULSES 


Irradiation  Apparatus 

The  irradiation  apparatus  used  in  these  studies  had  two  somewhat  different 
configurations:  the  first  designed  for  experiments  using  the  entire  mode- 
locked  pulse  train  of  ^100  pulses,  and  the  second  for  work  using  a  single 
pulse  selected  from  near  the  beginning  of  the  pulse  train. 

Configuration  Used  for  Pulse  Train  Studies--A  schematic  of  the  apparatus 
is  shown  in  Figure  4.  The  infrared  pulse  train  was  frequency-doubled  to  530  nm 
in  a  Type  I  phase-matched  KD*P  crystal  angle-tuned  for  maximum  conversion  ef¬ 
ficiency  (n.12%).  The  530-nm  component  was  separated  from  the  1060-nm  light  by 
a  dichroic  beamsplitter  (DBS).  A  Schott  KG-3  filter  blocked  residual  1R  from 
the  green  beam.  A  half-wave  retardation  plate  (a/2)  rotated  the  polarization 
of  2nd  harmonic  light  from  the  vertical  to  the  horizontal  plane.  An  uncoated 
pellicle  beamsplitter  (PL^)  sampled  a  small  fraction  (^1.7Z)  of  the  green  beam 
for  pulse  chronometry.  The  green  pulse  train  was  attenuated  to  approximately 
the  desired  energy  level  by  a  neutral  density  filter  stack  (NDF^).  An  Iris 
diaphragm  (ID)  selected  the  central  3  mm  from  the  green  beam,  which  had  di¬ 
verged  to  a  diameter  of  %10  mm  at  the  location  of  this  aperture. 

The  2nd  harmonic  pulse  energy  was  monitored  by  a  Laser  Precision  RkP-331 
pyroelectric  energy  probe  calibrated  by  the  manufacturer  using  standards  trace¬ 
able  to  the  National  Bureau  of  Standards.  The  pulse  energy  was  sampled  by  an 
uncoated  fused-silica  beamsplitter  (BS,  total  reflection  coefficient  =  7%, 
both  faces),  and  focussed  by  a  15-cm-f.l.  fused-silica  lens  into  the  aperture 
of  the  probe.  The  energy  in  the  pulse  train  was  displayed  on  a  Laser  Precision 
RkP-3230  digital  readout  unit. 

Additional  beam  attenuation  was  provided,  when  desired,  by  a  set  of  cali¬ 
brated  fused-silica  neutral  density  filters  (NDF^). 

Configuration  Used  for  Single-Pulse  Studies--The  setup  used  for  single-pulse 
irradiation  was  somewhat  different  from  that  described  above.  (See  Fig.  5.) 
First,  the  Pockels  cell  pulse-selecting  system,  installed  ahead  of  the  2nd 
harmonic  generator,  was  activated  as  described  in  Appendix  A.  Because  of  the 
low  energy  per  pulse  (^1%  of  the  total  energy  in  the  pulse  train),  the  half¬ 
wave  retardation  plate  (a/2)  was  placed  downbeam  from  pellicle  PL^.  This 
increased  the  reflection  coefficient  of  the  pellicle  by  about  an  order  of 
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Tigure  4.  Schematic  of  apparatus  for  irradiation  of  primate  eyes 
with  entire  trains  of  ultrashort  2nd  harmonic  (530  nm) 
pulses  derived  from  a  mode-locked  Nd:Glass  laser. 
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Figure  5.  Schematic  of  apparatus  for  irradiation  of  primate  eyes 
with  single  ultrashort  2nd  harmonic  (530  nm)  pulses 
derived  from  a  mode-locked  Nd:Glass  laser. 
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.tMyni  i-i>c]«  ..  i  iil>;  the  incident  pulse  was  now  vertically  polarized.  Without  this 
modification  not  enough  light  would  have  entered  the  pulse  chronometer  system 
to  yield  a  detectable  signal. 

The  most  significant  change  made  for  the  single-pulse  studies  was  the 
audition  of  a  beam-reducing  Galilean  telescope,  comprised  of  lenses  and  L-,. 
The  diameter  of  the  beam  was  reduced  by  a  factor  of  3.7  so  that  all  of  the 
green  pulse  energy  passed  through  the  3-mrn-aperture  iris  diaphragm.  This  beam 
reduction  was  necessary  so  that  the  single  pulses  could  reliably  trigger  the 
pyroelectric  energy  probe.  Without  the  beam  reduction,  only  about  10%  of  the 
pulse  energy  was  transmitted  by  the  aperture,  resulting  in  too  low  a  pulse 
energy  at  the  energy  probe  for  proper  triggering.  Unfortunately,  the  use  of 
the  beam-reducing  telescope  increased  the  beam  divergence  by  a  factor  equal  to 
the  beam  reduction  ratio.  Thus,  the  beam  divergence,  and  hence  the  retinal 
irradiation  spot  size,  for  the  single-pulse  experiments  was  3.7  times  greater 
than  in  the  experiments  utilizing  the  entire  pulse  train.  (The  method  used  to 
determine  the  beam  divergence  is  discussed  in  the  section  below.)  The  prob¬ 
lem  of  poor  energy  probe  triggering  on  single  pulses  was  not  discovered  until 
after  all  of  the  pulse  train  data  had  been  obtained,  so  the  two  sets  of  experi¬ 
ments  could  not  be  carried  out  using  the  same  beam  divergence. 

Control  of  pulse  energy  incident  on  the  eye  for  the  single-pulse  case 
was  effected  by  placing  calibrated  neutral  density  filters  (NDF^)  between  the 
energy-sampling  beamsplitter  and  the  eye  to  be  irradiated.  The  pulse  energy 
was  monitored  as  described  earlier. 

Apparatus  Common  to  Both  Confi gurations--The  pulse  chronometer  system, 
based  on  a  transverse-gated  optical  Kerr  effect  shutter,  is  identical  to  that 
used  in  our  macromolecul ar  irradiation  studies  (Appendix  A).  It  provides  an 
on-line  measure  of  the  pulse  duration  of  the  IR  laser  output.  The  frequency- 
doubled  pulse  duration  can  be  deduced  from  the  fact  that  its  intensity  scales 
as  the  square  of  IR  intensity.  Thus,  assuming  that  the  laser  pulses  have  a 
Gaussian  temporal  shape,  the  duration  of  the  530-nm  pulse  is  1//2  that  of  the 
IR  pulse.  In  this  work,  the  IR  pulse  duration  was  in  the  range  of  4-13  psec, 
with  an  average  of  9  psec.  It  follows  that  the  530-nm-pul se  durations  fell  in 
approximately  the  3-9-psec  range,  with  an  average  of  ^6  psec. 

The  beam  spatial  profiles  were  measured  by  directing  the  530-nm  beam  onto 
a  lensless  Fairchild  CCD  camera  having  1024  channels/inch.  The  linear  CCD 


diode  array  was  placed  at  the  same  distance  from  the  iris-diaphragm  aperture 
stop  as  was  the  pupil  of  the  test  primate.  The  beam  profile  was  approximately 

Gaussian  in  both  the  single-pulse  and  entire  pulse-train  configurations  with 

2 

a  1/e  diameter  of  <4  mm  at  the  pupil  in  each  case.  Beam  divergence  was  mea¬ 
sured  in  a  similar  manner,  with  the  CCD  array  placed  at  the  focal  point  of  a 
50-cm-f.l.  lens.  The  divergence  was  n-6  mrad  for  the  single  pulse  and  '>4.6  mrad 
for  the  pulse  trains,  indicating  irradiated  spot  diameters  on  the  retina  of 
''■78  pm  and  %2l  pm,  respectively,  based  on  a  typical  value  of  1.3  cm  for  the 
focal  length  of  the  primate  eye  (53). 

The  experimental  primates  were  mounted  in  a  prone  position  on  an  adjust¬ 
able  platform  having  3  degrees  of  translational  freedom  (x,y,z)  and  2  degrees 
of  angular  freedom  (azimuth  and  elevation).  The  pupil  of  the  subject  eye  could 
be  positioned  with  respect  to  the  laser  beam  axis  by  one  or  more  of  these 
mechanical  adjustments.  Visual  examination  of  each  primate  retina  was  per¬ 
formed  for  proper  placement  of  the  test  exposures  and  for  ore-  end  postexposure 
examination  to  determine  incurred  damage. 

A  dielectric-coated  pellicle  beamsplitter  ( PL^)  ■.  attached  to  a  two-posi¬ 
tion  swingaway  mount  affixed  to  the  fundus  camera  objective  barrel,  directed 
the  laser  pulse  into  the  eye.  In  the  "down"  position  the  beamsplitter  reflect¬ 
ed  21  0.5&*  of  the  incident  laser  energy  into  the  eye,  collinearly  with  the 
optic  axis  of  the  fundus  camera,  thus  permitting  simultaneous  observation  of 
the  irradiation  event.  This  on-line  observation  capability  permitted  precise 
positioning  of  the  irradiated  spot,  allowing  corrections  for  eye  motion  right 
up  to  the  time  the  laser  was  fired.  With  the  beamsplitter  in  the  "up"  posi¬ 
tion,  the  eye  could  be  examined  and  photographed  with  no  loss  of  illumination, 
or  vignetting. 

Experimental  Protocol 

Ten  Macaca  fasci cul ari s  primates  were  used  in  the  experiments.  Prior 
to  the  irradiation  studies  the  primates  were  refracted  in  each  eye  to 

*The  45*  reflection  coefficient  of  the  dielectric-coated  beamsplitter 
for  horizontally  polarized  light  at  530  nm  was  determined  experimentally, 
using  a  collimated  CW  light  source  of  the  correct  wavelength  and  polarization, 
chopped  at  a  frequency  of  330  Hz,  and  a  photomultiplier  detector.  The  re¬ 
flectivity  was  deduced  from  measurements  of  the  transmission  factor,  using 
the  fact  that  absorption  losses  in  the  coating  are  less  than  0.2%. 


do  tons i no  the  correct  i  ve  lenses  required.  (I’upi  1  dilition  and  cycioplegia  were 
obtained  by  instillation  of  two  drops  of  Kupfer's  solution,  a  one-to-one  mix¬ 
ture  of  1  cyclopentol ate  and  10  phenylephrine.)  Any  eye  with  a  refractive 
error  greater  than  1.50  diopters  in  any  meridian  was  not  used  in  the  experi¬ 
ments  . 

One  day  prior  to  the  retinal  irradiations  of  a  given  primate,  pupillary 
dilation  was  initiated  by  instillation  of  two  drops  of  atropine  sulfate  (4% 
solution).  This  was  followed,  immediately  prior  to  irradiation,  by  two  drops 
of  Kupfer's  solution.  Approximately  30  minutes  prior  to  irradiation,  each 
animal  was  deeply  sedated  with  an  i ntramuscul ar  injection  of  0.6  cc  ketamine 
HC1 .  Booster  doses  of  0.2  cc  ketamine  were  administered  during  the  course 
of  the  experiment  as  needed.  No  additional  general  anesthetics  or  tranquiliz¬ 
ers  were  required.  Eye  movements  were  eliminated  by  retrobulbar  injections 
of  xylocaine  (0.35  cc  each  side  of  the  orbit).  A  subcutaneous  injection  of 
0.15  cc  atropine  sulfate  served  to  suppress  drooling. 

During  the  exposure  procedures  the  eyelids  were  kept  retracted  by  a  pedi¬ 
atric  stainless-steel  speculum.  Corneal  dessication  was  prevented  by  frequent 
irrigation  with  normal  saline,  using  a  modified  hypodermic  syringe  equipped 
with  a  13-6  short  cannula  aimed  at  the  eye.  The  syringe  was  connected  to  the 
saline  solution  container  through  a  T-type  double  ball  valve.  The  syringe 
acted  as  a  positive  displacement  manual  pump  in  this  configuration  and  was 
adjusted  to  deliver  2  cc  of  saline  per  stroke  of  the  spring-1 oaded  plunger. 

To  guide  the  placement  of  the  irradiation  sites,  eight  marker  lesions, 
four  arranged  vertically  and  four  horizontally,  were  placed  adjacent  to  the 
macula,  as  shown  in  Figure  6.  These  markers  defined  a  cartesian  coordinate 
system  for  the  16  test  exposures  within  the  macula.  (In  some  cases  a  fifth 
row  or  column  was  also  irradiated  for  a  total  of  20  exposure  sites  within  the 
macula.)  The  horizontal  marker  row  was  always  inferior  to  the  macula,  whereas 
the  vertical  row  was  temporal  in  the  right  eye  and  nasal  in  the  left  eye.  The 
marker  lesions  were  produced  with  the  same  laser  system  used  to  produce  the 
experimental  lesions,  but  at  considerably  higher  energy.  The  entire  pulse 
train  of  530-nm  light  was  used  at  an  incident  energy  of  30-50  yJ.  The  result¬ 
ing  lesions  appeared  immediately  in  most  cases  and  had  a  diameter  ^3  times  the 
spot,  size  of  the  laser  beam  at  the  retina.  Most  of  the  marker  lesions  appeared 
a.  whitish  discolorations ,  although  a  few  exhibited  some  subretinal  hemorrhag¬ 
ing,  probably  a  result  of  damaging  hidden  capillaries. 
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Cross-hairs  in  the  eyepiece  of  the  fundus  camera  were  used  to  align  the 
exposure  site  with  the  marker  lesions.  The  animal  was  moved  relative  to  the 
laser  beam  to  change  the  exposure  site.  The  approximate  energy  range  of  a 
given  series  of  exposures  was  adjusted  by  the  first  neutral  density  filter 
stack  (NDF  ) .  (See  Fig.  4.)  Final  adjustments  were  made  with  the  calibrated 
neutral  density  filters  (NDF^).  The  laser  itself  provided  a  degree  of  random¬ 
ness  in  the  incident  energy  as  a  result  of  its  fluctuations  in  output  from 
shot  to  shot.  The  IR  output  of  the  laser  typically  fluctuates  over  a  range 
of  -25%  of  its  average  output.  Since  the  frequency-doubled  energy  scales  as 
the  square  of  the  IR  output,  the  530-nm  incident  energy  fluctuated  over  a 
range  of  50%  from  shot  to  shot.  In  cases  where  a  single  pulse  was  selected 
from  the  pulse  train,  the  fluctuation  was  even  greater  because  of  the  varying 
efficiency  of  pulse  selection  by  the  Pockels  cell  from  shot  to  shot.  (See 
Appendix  A.)  This  randomness  in  the  energy  delivered  to  successive  exposure 
sites  was  useful  in  preventing  bias  when  evaluating  the  postexposure  sites. 

Results 

The  maculae  of  each  primate  were  examined  with  the  fundus  camera  by  two 
observers  for  the  presence  or  absence  of  visible  lesions  at  1  hour  and  24  hours 
post  exposure.  A  lesion  was  defined  as  the  smallest  observable  circular  dis¬ 
coloration  (usually  whitish  or  light  gray)  differing  from  the  retinal  back¬ 
ground.  Using  this  criterion,  148  exposures  were  made  in  the  range  of  0.089- 
19.5  uJ  with  the  whole  pulse  train,  and  158  exposures  in  the  range  of  0.01- 
7.1  cJ  with  single  pulses,  to  establish  the  ED^q  point  (the  energy  required  to 
produce  a  lesion  in  50  of  the  events)  in  each  case. 

Before  the  data  were  analyzed,  the  pulse-train  and  pulse-width  data  were 
carefully  examined  to  eliminate  invalid  exposure  sites.  The  criteria  for  this 
test  were  as  fol lows: 

(1)  For  pulse  trains,  the  train  had  to  be  clean;  i.e.,  consist  of  a 
series  of  '100  equally  spaced  pulses.  Interlaced  multiple-pulse  trains  or 
the  presence  of  spurious  pulses  between  the  main  pulses  were  cause  for  data 
rejection. 

(2)  For  single  pulses,  more  than  90  of  the  energy  had  to  be  in  one 
pulse.  On  a  number  of  occasions  two  successive  pulses  in  the  train  were  switched 
out  by  the  Pockels  cel!.  Such  events  were  not  used  in  the  data  reduction. 
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Occasionally,  spurious  pulses  spaced  a  few  tens  of  picoseconds  from  the  select¬ 
ed  pulse  would  appear.  Those  data  shots  were  also  rejected. 

In  addition,  exposure  data  were  rejected  if  during  the  actual  laser  shot 
the  laser  spot  geometry  deviated  appreciably  from  a  circular  geometry  or  if  it 
appeared  fuzzy  and  significantly  larger  than  normal.  (Such  effects  could  be 
caused  by  gross  refractive  errors  resulting  from  distortion  of  the  cornea  by 
the  speculum  or  by  side  effects  of  the  retrobulbar  injections.)  Using  this 
criterion  alone,  all  data  from  three  eyes  had  to  be  rejected. 

In  the  case  of  the  whole  pulse  train  experiments,  77  exposure  events  passed 
the  above  criteria.  With  the  single-pulse  experiments,  108  events  passed. 
These  numbers  were  sufficient  to  establish  reliable  ED^q  thresholds  for  each 
category. 

The  data  were  submitted  to  standard  statistical  probit  analyses  (54)  to 
determine  the  ED<-q  points  (55).  The  analyses  were  carried  out  on  all  exposure 
sites  in  all  eyes  exposed  to  a  given  set  of  conditions.  Implicit  in  this  com- 
bined-probit  approach  is  the  assumption  that  the  variability  from  eye  to  eye  is 
no  greater  than  the  variability  among  sites  within  a  given  macula.  This  assump¬ 
tion  has  been  demonstrated  by  others  to  be  valid  (23).  The  calculations,  how¬ 
ever,  were  carried  out  separately  for  the  lesion/no-lesion  data  of  each  of  the 
two  observers.  Table  9  presents  the  1-hour  and  24-hour  ED5Q  points  for  each 
observer,  with  the  associated  95%  confidence  ranges  noted  in  parentheses. 


TABLE  9. 

RETINAL  DAMAGE 

THRESHOLDS  AT 

530  nm 

Energy  density  at  retina* 

Pul  se 

Si ngle 

Evaluation 

(Average  of  both  observers) 

p 

trai  n 

pul  se 

time 

Obs.  A 

Obs.  B 

( J/crn  ) 

X 

1  hr 

4.22  nJ 

4.46  yJ 

1.1 

(3.64-4.91) 

(3.76-5.29) 

X 

24  hr 

1.89  yJ 

2.31  y  J 

0.54 

(1.49-2.38) 

(1.87-2.84) 

X 

1  hr 

0.6!  J 

1.03  yJ 

1.5x10"^ 

(0.38-0.98) 

(0.65-1.64) 

X 

24  hr 

0.24  yJ 

0.24  ,:J 

4.4x10'° 

(0.17-0.35) 

(0.15-0.39) 

These  data 

corrected 

for  transmission 

factor  of  0. 

88  of  clear  ocular 

r  ed  i  a 

at  530  nm 

(29). 
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The  EDj-y  thresholds  for  single-pulse  lesions  are  considerably  lower  than 
the  thresholds  obtained  with  the  entire  pulse  train.  This  is  especially  evi¬ 
dent  in  the  24-hour  results  in  which  there  is  apparently  a  one-order-of-magni - 
tude  difference  between  the  single-pulse  and  whole  pulse-train  thresholds.  The 
difference  is  actually  much  greater  when  we  consider  the  fact  that  the  retinal 
spot  size  for  the  single-pulse  case  (^78  ;Jm)  is  m3. 7  times  greater  than  for  the 
pulse-train  case  (Ml  urn) .  Thus  the  energy  density  at  the  retina  is  an  addi¬ 
tional  factor  of  13.7  lower  for  the  single-pulse  case.  In  other  words,  the 
actual  threshold  for  damage  for  the  single-pulse  case  is  about  two  orders  of 
magnitude  below  that  for  an  entire  pulse  train,  in  terms  of  the  energy  density 
deposited  at  the  retina,  as  is  evident  in  the  last  column  of  Table  9. 

An  interesting  deduction  that  can  be  made  based  on  our  results  is  that, 
for  the  case  of  the  pulse  train,  the  lesion  event  may  actually  be  caused  by 
the  first  pulse  that  reaches  the  threshold  energy  for  a  single  pulse.  General¬ 
ly  this  occurs  near  the  beginning  of  the  train.  The  reasoning  is  simple: 

There  are  MOO  pulses  in  the  train  and  thus  the  energy  per  pulse  is  M/100  of 
the  total,  which  is  approximately  the  same  factor  by  which  the  single-pulse 
threshold  is  smaller  than  the  pulse-train  threshold  in  terms  of  energy  density 
at  the  retina.  The  remaining  pulses  in  the  train  thus  would  seem  to  only  add 
"insult  to  injury."  This  also  suggests  that,  if  no  pulse  in  the  train  meets 
the  threshold  energy  value  for  single  pulses,  there  may  be  no  damage  regardless 
of  the  number  of  pulses  in  the  train,  as  long  as  the  total  pulse-train  duration 
is  less  than  a  few  microseconds.  For  longer  trains,  thermal  effects  may  begin 
to  come  into  play  and  this  hypothesis  would  be  invalid. 

It  is  instructive  to  compare  the  retinal  damage  thresholds  at  530  nm  to 
the  results  of  our  earlier  vesicle  damage  studies  at  1060  nm,  in  terms  of  inci¬ 
dent  energy  density,  power  density,  and  corresponding  electric-field  strength. 
Table  10  shows  that  the  24-hr  retinal  thresholds  are  remarkably  close  to  the 
vesicle  thresholds  for  both  the  pulse  train  and  single  pulses,  also  that  the 
damage  threshold  for  a  single  pulse  is  about  two  orders  of  magnitude  below 
the'  for  i'i  enfin  pul  .<•  train  of  100  pulses  for  both  the  vesicles  and  the 
r<  r  iii  i .  :  lii  .  i,  ,t  to  ’urj  indication  that  the  damage  mechanism  in  both  cases 

i  de;i‘  r,  i»  nf  .n  :hc  .  tk  power  density,  and  hence  the  peak  electric  field  in 
trio  individual  :  ■  •  j  ■  ,  '  other  than  on  the  total  deposited  energy  density.  The 

1 1  a*  •••suit  ••deuce  to  the  hypothesis  forwarded  earlier  that  the 


first  pulse  in  the  train  to  reach  or  exceed  the  single-pulse  threshold  will 
cause  the  damage.  Subsequent  pulses  in  the  train  would  increase  the  number 
of  damaged  vesicles  or  cells  but  would  not  be  necessary  for  a  detectable  chang 
or  lesion. 

TABLE  10.  COMPARISON  OF  THRESHOLD  DAMAGE  DATA  FOR  VESICLES  AND  RETINAS 

Vesicle  Damage  Thresholds  Retinal  Damage  Thresholds* 

(  =  1060  nrn,  At  10  psec)  (‘  =  530  nm,  At  6  psec) 

Peak  elec- 
Energy  Power  trie  field 

J/cn GW/cm2  V/cm  x  ]QJ 

Pulse  0.6  0.6  5.8 

trai  n 

Single  9xl0"3  0.9  7.1  1  hr  1.5xl0~2  2.5  13 

pulse  24  hr  4.4X10*"5  0.73  6.4 

*Retinal  damage  thresholds  are  corrected  for  the  transmission  factor  of 
the  clear  ocular  media  at  530  nm.  T&30  =  0.88  (29). 

Since  the  electric  fields  associated  with  the  threshold  values  are  nearly 
an  order  of  magnitude  greater  than  the  membrane  potentials  of  both  the  vesicle 
and  cell  lipid  bi layer  membranes,  it  appears  reasonable  at  this  time  to  ascribe 
the  ocular  damage  to  membrane  disruption  by  electrostricti ve  forces. 

This  tentative  identification  of  the  damage  mechanism  must  be  tempered 
somewhat,  at  least  at  present,  by  a  number  of  factors.  First,  the  wavelengths 
used  in  the  vesicle  and  retinal  work  are  not  the  same.  Thus,  photobiologi cal 
processes  may  be  present  in  the  retinal  case  and  the  damage  may  not  be  a  result 
of  membrane  disruption  by  electric  fields.  Unfortunately ,  it  was  not  possible 
to  generate  sufficient  power  in  the  2nd  harmonic  of  the  Nd:Glass  laser  to 
attain  the  1060-nm  threshold  levels  for  the  vesicles  at  530  nm.  Such  an  expc  - 
ment  would  require  a  laser  amplifier  of  gain  3,  which  was  not  available  to  us 
at  the  time  of  the  research.  Thus,  the  wavelength  dependence  of  vesicle  Jan-jv. 
could  not  be  determined. 

Second,  the  vesicle  data  are  based  on  relatively  few  individual  experi¬ 
ments,  so  their  damage  threshold  determination  is  statistically  less  reliable 
than  for  the  retinae.  The  reason  for  the  relatively  few  data  points  is  that 
the  preparation  of  the  vesicle  samples,  their  irradiation,  and  their  sub- 
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Peak 
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fiel  d 
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GW/ cm2 

V/cm 

x  10b 

1.1 

1.8 

10 

0.54 

0.90 

7.1 

exam' nation  by  dynamic  light  scattering  are  very  time  consuming.  No 
more  than  one  or  two  samples  a  week  could  be  studied  in  this  manner,  whereas 
all  of  the  retinal  data  (320  data  points)  were  gathered  in  a  matter  of  a  few 
weeks . 

Finally,  the  vesicles  we  studied  were  not  subject  to  osmotic  pressure 
imbalances  or  chemical  potential  gradients  across  the  membrane,  whereas  such 
effects  are  present  in  living  cells.  In  the  living  cell  the  disruption  of  the 
membrane  by  the  electric  field  in  a  single  threshold  laser  pulse  would  be 
facilitated  by  the  osmotic  pressure  imbalance  across  the  membrane,  resulting 
in  irreversible  damage. 


Comparison  with  Other  Work 

The  24-hr  postexposure  results  reported  here  are  compared  in  Table  11 
with  the  results  of  similar  experiments  carried  out  by  Ham  et  al .  at  1064  nm 
(28),  Goldman  et  al.  at  both  532  nm  and  1064  nm  (29),  and  Taboada  et  al.  at 
1060  nm  (30,31).  A  considerable  discrepancy  appears  not  only  in  the  ED^g 
values  but  also  in  the  threshold  power  density  and  electric  field  at  the 
retina.  In  addition  to  differences  in  wavelength  and  pulse  duration,  a  number 
of  possible  factors  related  to  experimental  protocol  and  the  definition  of  a 
threshold  lesion  can  significantly  affect  the  experimental  results. 

One  factor  may  be  that  in  none  of  the  work  is  the  diameter  of  the  irradi¬ 
ated  spot  size  at  the  retina  known  with  a  great  degree  of  certainty.  This  is 
a  difficult  parameter  to  measure  in  vivo  (the  diameter  of  the  lesion  is  not 
necessarily  the  same  as  that  of  the  irradiated  spot);  consequently,  the  spot 
size  is  usually  estimated  from  a  knowledge  of  the  laser  beam  divergence  and 
an  estimate  of  the  effective  focal  length  of  the  emmetropic  primate  eye. 

Although  all  researchers,  including  ourselves,  have  used  corrective 
lenses  if  the  refractive  error  of  the  primate  eye  exceeded  j.0.5  diopter  in  any 
plane,  nonintrinsic  refractive  errors  can  arise  during  the  irradiation  experi¬ 
ments.  On  a  number  of  occasions,  for  example,  we  have  observed  otherwise 
normal  eyes  become  highly  astiqmatic  durinq  an  experiment,  possibly  as  a  re¬ 
sult  of  delayed  adverse  reaction  to  the  retrobulbar  injections  or,  to  a  lesser 
extent,  in  response  to  the  forces  exerted  by  the  eyelid  retracting  speculum. 

On  other  occasions  we  also  noted  that,  the  visibility  of  the  retina 
through  the  fundus  camera  was  greatly  reduced  due  to  opac i f ication  of  the 
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TABLE  11.  COMPARISON  OF  PRESENT  RETINAL  THRESHOLDS  WITH  RESULTS 
OF  OTHER  INVESTIGATIONS  (24  HR  POST  EXPOSURE) 


cornea.  Frequent  irrigation  with  normal  saline  hid  not  always  improve 
visibili ty. 

The  above  factors  can  lead  to  considerable  error  in  the  estimation  of 
the  irradiated  spot  size  and  to  significant  attenuation  of  the  incident  beam 
due  to  corneal  scattering.  In  either  case,  the  radiant  exposure  at  the  retina 
could  be  significantly  lower  than  expected,  resulting  in  high  values  of  FD^. 

The  experiments  reported  here  are  the  only  ones  in  the  picosecond  regime 
in  which  the  fundus  was  observed  during  actual  firing  of  the  laser.  The 
laser  was  fired  only  under  conditions  of  maximum  clarity  of  the  ocular  media. 

If  any  gross  astigmatism  was  evident  at  any  time  during  an  experiment,  the 
data  for  that  eye  were  rejected.  Not  all  previous  work  has  followed  such  pre¬ 
cautions. 

Other  variances  in  exposure  protocol  and,  particularly,  variances  in 
threshold  determination  can  also  significantly  affect  the  experimental  results. 
In  two  cases  (28,29)  the  definition  of  a  threshold  lesion  was  one  which  became 
just  visible  fundoscopically  within  24  hr  of  exposure.  No  statistical  analyses 
of  the  data  were  undertaken  by  those  workers,  so  the  determination  of  what  con¬ 
stituted  a  threshold  lesion  was  very  sensitive  to  the  observer's  subjective 
interpretation.  On  the  other  hand,  the  work  of  Taboada  et  al .  (30,31)  and  the 
results  presented  here  made  use  of  probit  analysis  to  determine  the  damage 
thresholds  and  are  thus  less  prone  to  subjective  errors. 

CORNEAL  DAMAGE  THRESHOLDS  INDUCED  BY  PICOSECOND  265-nm  LIGHT  PULSES 

Irradiation  Apparatus 

The  facility  used  in  the  UV  corneal  irradiation  experiments  was  in  most 
respects  similar  to  that  used  in  the  530-nm  retinal  exposure  experiments.  A 
schematic  of  the  apparatus  is  shown  in  Figure  7.  When  compared  with  the  facil¬ 
ity  shown  in  Figure  4,  the  UV  facility  differed  in  two  respects:  (1)  a  temper¬ 
ature  phase-matched  KD*P  crystal  was  added  to  frequency-double  the  53C-nm  light 
to  265  nm;  (2)  the  primate  was  placed  so  that  the  UV  beam  impinged  on  the 
cornea  directly,  without  deflection  by  a  beamsplitter  or  mirror. 

In  addition,  a  pair  of  Schott.  UG-5  UV  bandpass  filters,  having  a  combined 
transmission  of  50  at  265  nm  and  <10~^  at  530  nm,  was  placed  at  the  output  of 
the  4th  harmonic  crystal  to  eliminate  the  530-nm  light.  A  20-cm-f.l.  fused 
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quartz  (Suprasil)  lens  (L^)  was  used  to  vary  the  bean  diameter  impinging  on 
the  cornea  to  vary  the  incident  intensity  over  the  desired  ranqe.  A  set  of 
calibrated  UV  neutral  density  filters  (NDF)  was  used  for  additional  control 
of  intensity.  All  neutral  filters  ahead  of  the  UV  generator  were  removed  to 
assure  maximum  conversion  efficiency.  In  this  same  regard,  the  beam  diameter 
of  the  530-nm  light  was  reduced  by  a  factor  of  3.7  by  lens  pair  to  in¬ 

crease  the  intensity  and  thus  the  UV  conversion  efficiency  in  the  FHG  crystal. 
The  aperture  stop  was  set  at  3  mm  as  in  the  case  of  the  530-nm  experiments. 

The  beam  diameter  at  the  cornea  was  set  by  varying  the  distance  between 
the  cornea  and  the  focusing  lens  (L^).  Maximum  diameter  was  obtained  by  re¬ 
moving  the  lens.  The  cornea-to-lens  separations  used  in  the  experiments  were 
10,  15,  and  20  cm.  Other  values  of  incident  energy  density  were  obtained  by 
means  of  the  calibrated  UV  neutral  density  filters.  The  actual  beam  spot  size 
at  the  cornea  was  determined  by  a  fluorescence  technioue.  At  the  various  loca¬ 
tions  of  the  cornea  relative  to  the  focusing  lens,  a  strip  of  Scotch  "Magic" 
transparent  tape  coated  with  an  ethanolic  solution  of  coumarin  dye  (Eastman 
X5419)  was  placed  perpendicular  to  the  UV  beam.  The  very  thin  coating  of  dye 
left  on  the  tape  where  the  alcohol  evaporated,  fluoresced  strongly  without  any 
blooming.  A  closed-circuit  TV  camera  was  focused  on  the  rear  surface  of  the 
tape  to  record  the  fluorescent  spot.  The  intensity  distribution  in  this 
spot  was  determined  by  the  same  video  analyzing  system  used  to  measure  pulse 
widths  with  the  optical  Kerr  shutter.  Care  was  taken  to  eliminate  any  stray 
visible  or  IR  light  and  to  assure  that  the  fluorescence  response  was  linear, 
the  latter  verified  by  means  of  UV  neutral  density  filters.  The  spot  inten¬ 
sity  profiles  determined  in  this  manner  closely  followed  the  Gaussian  profile 
2 

with  1/e  diameters  of  2.9,  1.6,  1.1,  and  0.76  mm,  corresponding  to  no  lens, 
and  distances  of  10,  15,  and  20  cm  from  the  focusing  lens.  These  beam  sizes 
were  used  to  calculate  the  incident  energy  densities  at  the  irradiated  sites. 

Pulse  duration  measurements  were  carried  out  as  described  in  Appendix  A. 
Since  the  4th  harmonic  generation  scales  closely  as  the  fourth  power  of  the 
laser  output,  the  UV  pulse  durations  are  half  the  pulse  durations  of  the  1060-nm 
fundamental.  Thus,  in  this  set  of  experiments,  the  average  duration  of  each 
pulse  in  the  UV  pulse  train  was  ^4.5  psec. 
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Experimental  Protocol 

The  same  ten  primates  (Macaca  fascicularis)  used  in  the  530-nm  retinal 

irradiation  experiments  were  used  in  the  UV  corneal  irradiation  experiments. 

Prior  to  exposure  each  primate  was  prepared  as  described  previously  for  tUe 

530-nm  exposures,  with  the  exception  that  retrobulbar  injections  were  nut 

necessary  here.  The  corneas  of  each  animal  were  examined  with  a  hand-held 

ophthalmoscope  prior  to  the  exposure  experiment  to  locate  any  preexisting 

corneal  lesions.  Only  one  eye  had  a  preexisting  lesion  and  its  character 

was  recognizably  different  from  laser-induced  les'ons.  Four  sites  were  exposed 

on  each  cornea,  as  shown  in  Figure  8.  Only  entire  pulse  trains  were  used  in 

this  study,  due  to  the  relatively  low  UV  energy  available;  the  sinqle-pulse 

energy  was  too  low  to  trigger  the  energy  probe  reliably.  A  total  of  80  sites 

2 

were  exposed  over  an  incident  energy  density  range  of  1.3  to  70  mJ/cm  ,  using 
the  method  described  in  the  Irradiation  Apparatus  section  to  vary  the  energy 
density. 

Results  and  Discussion 

The  corneas  were  carefully  examined  with  an  ophthalmoscope  immediately 
after  exposure,  1  hr  post  exposure,  and  24  hr  post  exposure.  In  addition,  at 
24  hr  post  exposure,  the  corneas  were  examined  with  a  Nikon  slit-lamp  bio¬ 
microscope.  This  examination  was  done  both  in  white  light  and  with  cobalt 
blue  light.  In  the  latter  case,  sodium  fluorescein  dye  was  introduced  into 
the  eye  by  the  standard  technique  of  inserting  a  dye-impregnated  paper  strip 
under  the  lower  eyelid  and  then  spreading  the  released  dye  over  the  cornea  by 
blinking  the  eyelids  manually.  The  sodium  fluorescein  adheres  preferentially 
to  the  corneal  lesions  and  glows  greenish  yellow  under  cobalt  blue  illumina¬ 
tion.  This  was  the  most  sensitive  test  for  the  existence  of  threshold  lesion:,. 
The  results  were  recorded  photographically  as  well  as  visually.  No  lesions 
were  detected  either  immediately  or  1  hr  post  exposure,  even  at  the  highest 
energy  densities.  However,  epithelial  lesions  were  observed  in  about  half 
the  cases  at  24  hr  post  exposure.  Under  white-light  illumination,  the  lesion- 
appeared  as  small  circular  grainy  opacifications,  having  a  diameter  somewhat 
smaller  than  the  UV  beam  diameter.  None  of  the  lesions  appeared  to  have  ary 
significant  depth.  Stained  with  sodium  fluorescein  and  observed  under  cobal; 
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blue  illumination,  the  lesions  stood  out  more  clearly.  In  all  cases  the 
lesion/no-lesion  events  were  recorded  by  two  observers. 

The  lenses  of  eyes  that  exhibited  corneal  lesions  were  also  examined  with 
the  slit-lamp  for  the  possible  presence  of  UV-induced  cataracts.  However, 
unlike  work  reported  at  longer  wavelengths,  >325  nm  (26,27),  no  signs  of  any 
lenticular  opacification  were  seen  in  any  of  the  eyes.  Retinal  examinations 
were  not  made  since  the  subject  animals  had  been  previously  used  in  the  530-nm 
retinal  damage  studies. 

The  data  were  statistically  analyzed  using  the  method  of  probits,  as  for 
the  530-nm  retinal  exposures.  An  ED^g  threshold  was  definable  only  24  hr  post 
exposure,  as  no  lesions  were  observed  1  hr  post  exposure.  Again,  only  data 
that  represented  well-defined  operation  of  the  laser  were  used  in  the  analysis. 
(A  total  of  61  events  were  usable.)  The  mean  of  observer  values  averaged  over 
the  eyes  gave  the  following  estimate  with  the  associated  95%  confidence  limits: 
EDrn  =  8.2  inJ/cm2  (6.3  -  10.7)  at  24  hr  post  exposure.  The  corresponding  peak 

bU  2 

power  per  pulse  in  the  train  at  the  cornea  was  18  MW/cm  ,  and  the  electric 
field  per  pulse  was  '.10s  V/cm. 

There  are  no  data  on  UV  laser  effects  to  which  these  results  can  be  com¬ 
pared  directly.  The  shortest  UV  laser  wavelength  previously  studied  for 

corneal  damage  effects  was  337  nm  from  a  pulsed  nitrogen  laser  producing  10-nsec 

2 

pulse  durations  (27).  In  that  work  the  damage  threshold  was  8.7  J/cm  ,  or 
three  orders  of  magnitude  higher  than  the  present  value.  No  correlations  can 
be  made  with  that  result  because  of  the  vastly  different  pulse  durations  and 
because  corneal  tissue,  as  all  biological  tissue,  exhibits  a  considerable 
variation  in  its  adsorption  spectrum  in  the  region  between  200  nm  and  400  nm 
(2,3). 

The  only  data  available  at  a  wavelength  close  to  the  265  nm  studied  hero 
are  those  of  Pitts  et  al.  (25)  at  270  nm  for  a  noncoherent  UV  source  (an  arc 
monochromator)  and  essentially  CW  exposure.  Pitts  and  his  coworkers  found 
that  the  peak  sensitivity  of  the  cornea  to  UV  photokeratitis  was  at  270  nm  for 
both  the  primate  and  human  eye,  and  they  deterurned  the  threshold  to  he 

O 

4  m,)/cm‘‘  for  both  subjects.  This  value  is  comparable  in  magnitude  to  our 
results  at  265  nm;  however,  except,  for  some  punctate  staining,  the  experi¬ 
ments  of  Pitts  did  not  produce  well-defined  corneal  opacities  similar  to  ours. 
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The  primary  damage  mechanism  for  picosecond  pulses  at.  265  nn;  is  somewhat 

speculative  at  this  time,  but  it  appears  to  be  photochemical  in  nature. 

Neither  the  power  density  nor  electric  field  per  pulse  in  the  mode-locked 

pulse  train  appears  to  be  sufficient  for  nonlinear  or  direct,  electric-field 

5 

effects.  For  example,  the  peak  electric  field  of  10  V/cm  carried  by  the 
individual  pulses  is  too  close  in  magnitude  to  the  electric  field  within  the 

c 

cell  membranes  (0.6x10  V/cm)  to  have  a  significant  effect. 

Although  the  mechanism  of  the  effect  of  UV  radiation  on  biological  sys¬ 
tems  is  only  generally  understood,  we  believe  that  the  most  likely  damage 
mechanisms  at  265  nm  are  the  photochemical  alteration,  such  as  denaturation 
and  coagulation,  of  proteins  and  nucleic  acids.  Especially  vulnerable  are 
unconjugated  nucleoproteins  of  the  cell  and  the  DNA  in  the  chromosomes.  Both 
are  particularly  susceptible,  with  maximum  sensitivity  at  265  nm,  since  this 
wavelength  lies  at  or  near  the  center  of  their  action  spectra  (2,3).  Further 
work  needs  to  be  carried  out  at  265  nm  in  a  manner  similar  to  our  experiments 
in  the  visible  and  near  IR  to  confirm  this  damage  mechanism. 
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APPENDIX  A 


PICOSECOND  LASER  IRRADIATION  FACILITY 

The  ultrashort-pulse  laser  facility  assembled  for  the  purpose  of  irradi¬ 
ating  selected  macromolecular  samples  and  primate  eyes  is  illustrated  in 
Figure  A-l.  It  consists  of  four  subsystems:  a  mode-locked  NdrGlass  laser, 
a  high-speed  Pockels  cell  pulse-switching  system,  a  pulse  chronometer  and 
video  detection  system,  and  a  pulse-energy  measurement  system.  Each  of  these 
is  described  below. 


Nd:Glass  Laser 

The  mode-locked  Nd:Glass  laser  consists  of  a  water-cooled,  Brewster-angled, 
Owens-Illinois  ED-2  glass  rod,  1.3 -cm  dia  x  22.9-cm  length,  pumped  by  two  EG&G 
linear  flashlamps  in  a  double  elliptical  reflector  cavity.  The  resonator  is 
formed  by  a  flat  99.7%  rear  reflector  (Mj)  and  a  10-m  radius  35%  output  reflec¬ 
tor  (M2).  Mode-locking  is  accomplished  by  a  flowing  ]-mm-thick  dye  cell  placed 
in  direct  contact  with  the  rear  reflector.  The  dye  solution  consists  of  East¬ 
man  9860  dye  in  dichloroethane  at  a  concentration  that  results  in  a  small - 
signal  transmission  factor  of  '^60/'  at  1060  nm.  An  iris  diaphragm  (IDj)  is  used 
to  control  transverse  mode  size  and  purity.  By  closing  it  down  to  an  aperture 
of  5-mm  dia  or  less,  TEMqo  output  can  be  obtained.  In  the  above  configuration 
the  laser  produces  a  train  of  'vlOO  horizontally  polarized  pulses  at  >  =  1060  nm, 
each  of  slO-psec  duration  and  >100-MW  peak  power,  spaced  at  5.6-nsec  intervals 
(the  round-trip  cavity  time). 

The  choice  of  dye  cell  geometry  is  critical  to  laser  performance.  We 
have  experimented  with  various  types  of  discrete  dye  cells  and  with  the  con¬ 
tacted  type,  with  and  without  dye  circulation,  and  have  found  that  the  circulat¬ 
ing  contacted  cell  produces  the  most  consistent  and  reproducible  pulse  trains, 
with  excellent  suppression  of  satellite  pulses.  Beam  stability  and  mode  purity 
are  also  optimized  with  this  configuration. 

At  one  point  in  the  program  we  experienced  difficulties  in  proper  mode¬ 
locking  of  our  glass  laser  on  account  of  the  poor  quality  laser-grade 
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I  ..I  dichloroethane  supplied  by  Eastman  Organic  Chemicals.  Apparently,  high 
quantities  residual  HC1  were  present  i n  their  ther.-avai  1  able  lot.  No.  A6B, 
which  destroyed  the  9860  saturable  dye  immediately  upon  mixing.  No  mode¬ 
locking  at  all  was  obtained  using  this  particular  lot  of  solvent.  A  search 
through  various  chemical  distributors  turned  up  some  dichloroethane  of  a  differ¬ 
ent  lot  number  (A4C),  which  worked  marginally.  The  dye  solution  mode-locked 
the  laser  successfully  for  two  or  three  consecutive  days,  but  then  it  degraded 
and  had  to  be  replaced. 

Because  of  these  difficulties,  for  some  of  the  experiments  the  glass  laser 
was  mode-locked  by  Eastman  9740  dye  in  chlorobenzene  instead  of  the  usual  9860 
dye  in  dichloroethane.  With  the  9740  dye  the  laser  produced  the  same  average 
total  energy  per  pulse  train,  but  the  number  of  pulses  in  a  train  was  at  most 
only  .30-50  compared  to  the  average  'x.SO-lOO  obtained  with  9860  dye.  Operation 
of  the  laser  was  much  more  erratic  with  9740  than  with  9860.  Multiple  and 
spurious  pulsing  occurred  on  at  least  half  of  all  mode-locked  shots.  Only  one 
in  three  or  four  firings  of  the  laser  resulted  in  mode-locking  at  all.  Further- 
niore,  due  to  the  higher  peak  power  obtained  with  9740  dye,  some  of  the  optical 
components  sustained  some  surface  or  bulk  damage.  Because  of  these  difficulties, 
the  use  of  9740  dye  was  discontinued. 

The  problem  of  impure  dichloroethane  was  finally  solved  by  filtering  the 
solvent  through  a  4"  column  of  basic  alumina  powder,  following  the  suggestion 
of  W.  Robinson  at  Texas  Tech  University.  The  purified  solvent  was  then  filtered 
free  of  residual  alumina  particles  by  means  of  a  0.45-um  Millipore  filter. 

The  entire  purification  process  is  slow  but  very  effective.  Excellent  mode¬ 
locking  has  been  achieved  using  the  purified  dichloroethane  as  the  solvent  for 
the  9860  dye.  Only  half  the  usual  concentration  of  dye  is  required,  and  the 
solution  remains  stable  for  weeks  at  a  time.  Mode-locking  reproducibility  has 
improved  considerably  also.  Spurious  or  multiple  pulsing  occurs  no  more  than 
once  in  20-30  shots.  Although  power  output  is  somewhat  lower  than  in  the  past 
(150-200  MW/cm  compared  to  200-300  MW/cm  ),  the  pulse  duration  is  less. 

Pulses  at  the  beginning  of  the  pulse  train  have  been  measured  to  have  durations 
in  the  6-10-psec  range,  whereas  previously  they  were  in  the  10-15-psec  range. 


l1ockels  Cell  Pulse-Switching  System 

To  permit  irradiating  the  macromolecular  samples  or  primate  eyes  with 
single  ultrashort  pulses  as  well  as  with  entire  pulse  trains,  a  provision  for 
switching  out  a  single  pulse  from  the  mode-locked  train  has  been  incorporated. 
A  high-speed  electro-optical  shutter  is  used,  consisting  of  two  crossed  thin- 
film  polarizer  pairs  (TFP)  on  either  side  of  a  high-speed  Lasermetrics  1071-FY 
Pockels  cell,  which  is  switched  by  a  ^.7-kV  pulse  provided  by  a  Lasermetrics 
type  8601  avalanche-transistor  Krytron-triggered  Blumlein  pulser.  The  Pockels 
cell  is  connected  to  a  50-n  terminator  via  a  30.5-m  (100-ft)  length  of  RG-8/U 

coaxial  cable.  The  thin-film  polarizers  are  at  Brewster's  angle  (56.5°)  and 
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are  stacked  in  pairs  to  yield  a  polarization  ratio  of  %2.8xl0  for  each  pair. 
The  output  of  the  laser  is  horizontally  polarized.  For  convenience  the 
polarization  vector  is  rotated  into  the  vertical  plane  by  a  half-wave  retarder 
plate  ( x/2 ) .  The  first  thin-film  polarizer  rejects  any  residual  horizontal 
polarization  component.  After  passing  through  the  inactive  Pockels  cell,  the 
first  pulse  in  the  train  is  totally  reflected  by  the  second  polarizer  stack  to 
an  ITT  FW  4014  biplanar  vacuum  photodiode,  whose  output  is  displayed  on  a 
Tektronix  type  519  oscilloscope. 

The  first  laser  pulse  triggers  the  oscilloscope.  Approximately  45  nsec 
later  a  fast-rising  step-function  voltage  pulse  appears  at  the  "+  Gate"  output 
of  the  oscilloscope.  This  signal  is  used  to  trigger  the  high-voltage  pulser 
which  activates  the  Pockels  cell.  The  +  Gate  output  delay  can  be  continuously 
varied  from  45  to  80  nsec,  thus  permitting  precise  timing  of  the  pulse  deliv¬ 
ered  to  the  Pockels  cell.  The  pulser  itself  has  an  additional  variable-delay 
control,  which  can  be  set  for  pulse  delays  of  ^100-400  nsec  relative  to  the 
triggering  signal  if  desired.  The  direct  mode,  which  affords  a  shorter  delay 
of  -v35  nsec,  has  been  used  in  our  experiments.  The  laser  pulse  that  happens 
to  pass  through  the  Pockels  cell  while  it  is  switched  on  has  its  polarization 
rotated  into  the  horizontal  plane  and  thus  passes  through  the  second  polarizer 
unimpeded.  The  remaining  pulses  in  the  train  arrive  in  the  Pockels  cell  after 
the  switching  pulse  and  hence  are  totally  reflected  by  the  second  polarizer. 

We  have  been  able  to  switch  out  clean  single  pulses  with  up  to  90"  effi¬ 
ciency,  with  a  good  degree  of  reproducibility.  Generally,  about  3  out  of  5 
shots  are  successful  (i.e.,  >70%  throughput,  with  no  measurable  bleed-through 


of  adjacent  pulses),  provided  the  mode-locked  pulse  train  is  devoid  of  any 
spurious  or  satellite  pulses  which  could  cause  premature  triggering. 

Pulse  Chronometer  System 

Temporal  width  measurement  of  the  selected  ultrashort  pulse  is  carried 
out  by  means  of  a  picosecond  streak  shutter  that  we  have  used  extensively  in 
past  work  (34,46,56,57).  The  selected  single  pulse  generates  2nd  harmonic 
light  (530  nm)  in  a  KDP  crystal  tuned  to  yield  an  SHG  conversion  efficiency 
of  about  1%.  The  superimposed  infrared  (IR)  and  green  pulses  are  separated 
at  the  dichroic  beamsplitter  (DBS)  (Fig.  A-l).  For  the  macromolecular  studies 
a  50%  beamsplitter  ( BS ^ )  directs  half  the  IR  pulse  energy  to  the  ultrafast 
streak  shutter.  (For  the  primate  experiments  all  the  IR  energy  is  directed 
into  the  streak  shutter,  since  the  irradiation  of  the  eye  is  carried  out  with 
the  2nd  or  4th  harmonics).  The  polarization  of  this  pulse  is  rotated  into  the 
vertical  plane  by  a  half-wave  retarder  (a/2).  Lenses  and  reduce  the 
beam  diameter  by  a  factor  of  four.  The  pulse  then  traverses  a  quartz  cell 
filled  with  carbon  disulfide  (CS^).  This  cell  is  located  between  two  high- 
quality  crossed  polarizers  (PI,  P2),  whose  polarization  axes  are  inclined  at 
45'  with  respect  to  the  polarization  of  the  IR  pulse.  These  three  components 
constitute  the  ultrafast  streak  shutter.  As  it  travels  through  the  CS^,  the 
IR  pulse  induces  a  narrow  zone  of  birefringence  in  its  immediate  vicinity  (46). 
To  an  observer  viewing  the  shutter  at  right  angles  to  the  IR  path,  the  effect 
is  that  of  a  narrow  "slit"  moving  across  the  line  of  sight  at  the  speed  of 
light  in  CS^  (1.84xl010  cm/sec).  The  shutter  thus  produces  a  streak  record 
of  light  pulses  incident  at  right  angles  to  the  IR  path. 

The  530-nm  pulse  split  off  at  the  dichroic  beamsplitter  is  directed 
toward  the  shutter  by  a  right-angle  prism  (PR)  and  expanded  horizontally  by  a 
pair  of  cylindrical  lenses  (L3,  L4j  to  illuminate  the  entire  length  of  the 
shutter,  where  it  is  sampled  by  the  IR  gating  pulse.  In  the  case  of  the  pri¬ 
mate  experiments,  only  a  portion  of  the  530-nm  light  is  directed  to  the  shut¬ 
ter,  by  means  of  a  pellicle  beamsplitter.  The  signal  exiting  from  the  shutter 
is  a  cross-correlation  between  the  gating  and  green  pulses  (57).  If  the  depth 
of  the  IR  gating  pulse  is  small,  the  transverse  dimension  of  the  transmitted 
green  pulse  is  essentially  the  same  as  the  geometrical  pulse  length  of  the  IR 
pul se  in  air. 


The  pulses  gated  by  the  shutter  are  detected  and  processed  by  the  video 
detection  and  display  system  (VDDS)  shown  in  Figure  A-2  (46,57).  The  shutter 
output  is  imaged  by  a  Telemation  TMC-1100  CCTV  camera  equipped  with  an  RCA 
4532A  silicon  vidicon  tube.  The  video  signal  is  processed  by  the  video  dis¬ 
play  control  unit  (VDC,  built  in-house)  and  displayed  on  an  RCA  CCTV  monitor. 
Superimposed  on  the  display  is  a  bright  rectangular  frame  generated  by  the 
VDC.  The  frame  height  can  be  varied  from  1  to  64  TV  lines,  and  its  width 
from  an  equivalent  of  64  lines  to  full  screen  width.  The  intensity  profiles 
of  the  TV  lines  within  the  frame  are  displayed  on  an  oscilloscope.  The  VDDS 
is  operated  in  the  single-shot  mode,  wherein  only  a  single  sweep  of  the  vidi¬ 
con  and  display  oscilloscope  occurs.  In  this  case  only  half  the  field  of  TV 
lines  is  swept;  i.e.,  only  the  odd-  or  even-numbered  lines.  This  avoids 
charge  leakage  from  the  transient  image  on  the  vidicon  in  the  time  between 
sweeps  of  the  odd  and  even  fields.  Thus,  in  this  mode  up  to  32  alternate 
lines  can  be  examined.  A  trigger  output  pulse  from  the  VDC  fires  the  laser 
at  the  start  of  the  vidicon  sweep. 

In  our  studies  we  have  used  a  frame  height  of  one  line  and  positioned  it 
to  provide  a  horizontal  slice  through  the  center  of  the  cross-correlation 
pulse  image.  The  oscilloscope  display  is  thus  a  plot  of  pulse  intensity  as  a 
function  of  line.  All  pulse  durations  quoted  are  measured  full  width  at 
half-maximum  and  are  deconvoluted  for  the  effect  of  the  finite  thickness  of 
the  birefringent  zone  in  the  shutter  medium. 

Pulse  Enerqy  Measurement 

The  IR  laser  pulse  energy  is  monitored  by  a  laser  Precision  RkP-331 
pyroelectric  energy  probe  and  RkP-3230  digital  display  unit.  This  system  was 
calibrated  by  the  manufacturer  using  standards  traceable  to  the  National 
Bureau  of  Standards.  A  20-cm-f.l.  lens  placed  10  cm  in  front  of  the  probe 
focuses  the  incident  IR  pulse  to  <3-mm  diameter  for  acceptance  by  the  probe 
aperture.  The  incident  pulse  is  sampled  by  means  of  an  uncoated  glass  beam¬ 
splitter  (BS^  in  Fig.  A-l). 
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Figure  A-2  .  Schematic  of  video  detection  and  display  system. 


APPENDIX  B 


EXPERIMENTAL  PROTOCOL  FOR  DNA  STUDIES 

Calf- thymus  DNA  solutions  were  carefully  prepared  at  a  1-mg/ml  con¬ 
centration  in  aqueous  sodium  chloride  (NaCl)  by  dissolving  the  DNA  in  a  cold 
room  (5°C)  for  periods  of  3-5  days,  at  a  carefully  controlled  stirring  rate 
of  ^1  cycle/sec.  Some  samples  included  the  Ca++  sequestering  agent,  EDTA; 
others  did  not.  In  each  case  both  the  test  and  control  samples  were  treated 
identically  except  for  irradiation  by  the  NdrGlass  laser.  Part  of  the  sample 
solution  was  used  to  fill  the  selected  test  cell  for  irradiation  by  the  pico¬ 
second  laser  pulses.  An  identical  cell  was  fill ea  at  the  same  time  to  act 
as  the  control  sample.  The  test  cell  was  filled  carefully  to  minimize  the 
occurrence  of  air  bubbles  at  the  optical  windows  or  in  the  optical  path  of 
the  laser  pulse. 

Two  types  of  UV-grade  fused-quartz  test  cells  were  manufactured  by  Preci¬ 
sion  Cells,  Inc.  for  the  irradiation  experiments.  One  cell  type  has  a  5 -mm 
ID  and  a  2-cm  path  length.  Two  filler  spouts  at  either  end  are  provided  for 
filling  and  flushing  the  cell.  It  is  referred  to  as  the  "short"  cell.  The 
other  type  is  similar  in  construction  to  the  short  cell  but  has  a  2-mm  ID 
and  a  10-cm  length  and  is  referred  to  as  the  "long"  cell.  It  is  used  when 
higher  incident  energy  levels  are  produced  by  focussing  down  the  incident 
beam  diameter.  In  each  case  the  incident  laser  pulse  just  fills  the  bore 
of  the  cell  and  thus  irradiates  the  entire  sample. 

Test  samples  were  irradiated  by  single  mode-locked  pulses,  entire  pulse 
trains,  and  successions  of  several  pulse  trains.  To  achieve  higher  energy 
densities,  the  smal 1 -diameter  long  cell  was  used.  The  test  and  control  samples 
were  then  taken  out  of  their  respective  cells  and  diluted  with  1-M  NaCl  solu¬ 
tion  in  the  ratio  of  0.5  ml  DMA  to  9.5  ml  NaCl.  The  solutions  were  filtered 
through  either  a  3.0-  or  1.25-um  filter  and  collected  in  scattering  cells 
for  the  dynamic  light  scattering  experiments.  Sample  handling  for  the  gel 
electrophoresis  and  low-shear  viscometry  is  discussed  in  the  main  text. 
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APPENDIX  C 


DYNAMIC  LIGHT  SCATTERING  FACILITY 

The  dynamic  light  scattering  apparatus  consists  of  a  CW  He-Ne  laser  oper¬ 
ating  at  632.8  nm  with  approximately  50-mW  power,  the  optical  detection  system, 
photon  counting  and  correlating  electronics,  and  a  PDP-12  computer  for  data 
analysis  and  storage.  The  block  diagram  of  the  experimental  setup  is  shown 
in  Figure  C-l.  The  laser  and  optical  detection  systems  are  mounted  on  a  vibra- 
tionally  damped  table  which  consistsofa  1360-kg,  3.7-m-long,  83-cm-wide- 
flange,  steel  I-beam  sitting  upon  16  free-floating  springs  grouped  into  four 
groups  of  four.  The  resonance  frequency  of  the  beam  with  respect  to  the  floor 
is  between  1  and  2  Hz,  as  per  design. 

The  output  of  the  laser  is  passed  through  an  optical  polarizer  and  a 
10-crn-f.l.  lens  to  focus  the  laser  to  its  minimum  waist  in  the  scattering 
cell.  The  optical  detection  system  consists  of  a  low-noise-selected  ITT-FW 
130  photomultiplier  tube  with  a  mu-metal  shield  enclosed  in  a  specially  con¬ 
structed  aluminum  housing  that  contains  the  photomultiplier  dynode  bias  elec¬ 
tronics.  The  optical  detection  system  is  mounted  on  a  triangular  optical 
rail  rigidly  attached  to  a  rotary  milling  table. 

A  portion  of  the  scattered  light  from  the  sample  passes  through  a  10-cm 
lens  placed  such  that  the  light  passing  through  the  lens  forms  a  divergent 
cone.  This  is  easily  achieved  by  placing  the  lens  so  that  the  distance  be¬ 
tween  the  focused  laser-beam  waist  in  the  sample  and  the  lens  is  less  than 
the  focal  length  of  the  lens.  A  series  of  apertures  between  the  photomulti¬ 
plier  and  the  lens  admit  light  from  only  a  rather  small  solid  angle.  The 
actual  collection  solid  angle  is  ultimately  determined  by  the  dimension  (0.25 
mm)  of  the  photoactive  cathode  and  the  divergence  angle  of  the  beam,  and  is 
slightly  less  than  the  coherence  solid  angle.  (Collecting  a  greater  solid 
angle  simply  includes  more  independently  fluctuating  K-vectors,  reducing  the 
apparent  signal-to-noise  ratio.) 

The  photomultiplier  output  is  fed  to  the  PAR  ampl i fier-discrimi nator , 
which  selects  the  photoelectron  pulses  and  transfers  them  to  the  pulse  invert¬ 
er.  This  in  turn  converts  them  to  TTL  pulses  suitable  for  input  to  the 
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nental  arrangement  for  dynamic  light  scattering. 


digital  photon-correlator.  The  output  correlation  function  is  stored  in  the 
memory  channels  of  the  Nicolet  signal  averager,  which  car  be  transferred  to 
the  PDP-12  computer  via  its  A-D  converter. 

The  digital  photon-correlator  is  a  256-channel  Chen-type  digital  clipped 
correlator  (DCC)  (58)  constructed  in-house.  A  block  diagram  of  the  DCC  is 
shown  in  Figure  C-2. 
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Figure  C-2.  Block  diagram  of  digital  clipped  correlator. 


